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FOREWORD 


The objective of the Tether Dynamics Simulation Workshop was to 
provide a forum for the discussion of the structure and status of 
existing computer programs which are used to simulate the dynamics of a 
variety of tether applications. A major topic of the workshop concerned 
itself with the purpose of having different simulation models and the 
process of validating them. Guidance on future work in these areas was 
obtained from a panel discussion consisting of resource and technical 
managers and dynamic analysts in the tether field. The conclusions of 
this panel are presented in Section HI, Workshop Summary and Panel 
Discussion Results. 

NASA Headquarters, Office of Space Flight, Advanced Programs, along 
with NASA Marshall Space Flight Center, sponsored the workshop. General 
Research Corporation and the American Institute of Aeronautics and Astro- 
nautics provided the workshop coordination. General Research Corporation 
also prepared the workshop proceedings. The workshop was held at the 
Hyatt Regency Crystal City Hotel in Arlington, Virginia in conjunction 
with the International Conference on Tethers in Space, which was held 
Sept. 17-19, 1986. 


PRECEDING PAGE BLANK NOT FU*» 


iii 


WORKSHOP ADMINISTRATIVE PERSONNEL 


WORKSHOP ORGANIZERS 

Charles C. Rupp William A. Baracat 

NASA Marshall Space Flight Center General Research Corporation 


ADMINISTRATIVE COORDINATORS 

Edward J. Brazill 
NASA Headquarters 


Terrence G. Reese 

General Research Corporation 


EXECUTIVE COMMITTEE 

Paul A* Penzo Peter M. Bainum 

Jet Propulsion Laboratory American Astronautical Society/ 

Howard University 

Charles C. Rupp 

NASA Marshall Space Flight Center William A. Baracat 

General Research Corporation 

Silvio Bergamaschi 
University of Padova, Italy 


PANEL MEMBERS 


Howard Flanders 

Martin Marietta Denver Aerospace 
Dave Arnold 

Smithsonian Astrophysical Observatory 
Joe Carroll 

Energy Science Laboratories 
Enrico Lorenzini 

Smithsonian Astrophysical Observatory 
Peter Bainum 

American Astronautical Society/ 

Howard University 


Henry Wolf 

Analytical Mechanics Associates 

Silvio Bergamaschi 
University of Padova, Italy 

Charles C. Rupp 

NASA Marshall Space Flight Center 
Paul Penzo 

Jet Propulsion Laboratory 

Bill Djinis 
NASA Headquarters 


WORKSHOP COORDINATION 

Pamela Edwards Mereille Gerard 

American Institute of American Institute of 

Aeronautics and Astronautics Aeronautics and Astronautics 


v 


PRECEDING PAGE BLANK NOT FILMED 



TABLE OF CONTENTS 


Page 

Foreword iii 

Workshop Administrative Personnel v 

I. WORKSHOP PRESENTATIONS 1 

"Tether Dynamics Simulation," A.K. Misra and V.J. Modi 3 

"Martin Marietta Simulation - HIFITS (Hi-Fidelity- 
Tethered-Satellite) Simulation," H. Flanders 25 

"The Tethered Satellite System on the Systems 

Engineering Simulator," R. Humble 57 

"GTOSS - Generalized Tethered Object Simulation 

System," D. Lang 73 

"Tethered Satellites Simulation," J.R. Glaese 95 

"Skyhook Program," D. Arnold 113 

"Slack3," G. Gullahorn 125 

"Artificial Gravity Laboratory," E. Lorenzini 141 

"Elasticity Effects in Tether Dynamics," S. Bergamaschi 161 

"TSS Subsatellite Attitude Dynamics and Control 

Laws Verification Programs," F. Venditti 187 

"Validation of TSS Engineering Simulation," Z. Galaboff 219 

"TSS-1 Dynamics Flight Experiments," G. Gullahorn 235 

"Tether Simulation Design for Mission Planning 

and Analysis," R. Deppisch and Y. Jani 249 

II. ADDITIONAL PRESENTATIONS 277 

"Preliminary Analysis of Damping Effect on Tethered 
Satellite," G. Bianchini 279 

"Rigorous Approaches to Tether Dynamics in Deployment 

and Retrieval," E. Antona 315 

III. WORKSHOP SUMMARY AND PANEL DISCUSSION RESULTS 331 

APPENDICES 

A. Workshop Agenda 337 

B. Workshop Attendees 341 


vii 


PRECEDING PAGE BLANK NOT FILMED 



I 

WORKSHOP PRESENTATIONS 


1 



Tether Dynamics Simulation 

Arun K. Misra 
McGill University 

Vinod J. Modi 

University of British Colombia 


RECEDING fAGE BLANK NOT FILMED 


3 
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PURPOSE OF SIMULATION 


I. TO STUDY THE GENERAL DYNAMICS OF THE SHUTTLE-TETHER-SUBSATELLITE 

SYSTEM DURING (i) STATIONKEEPING, 

(ii) DEPLOYMENT OF THE SUBSATELLITE, 

(m) RETRIEVAL OF THE SUBSATELLITE. 

II. TO DETERMINE APPROPRIATE CONTROL LAWS AND SIMULATE THE CONTROLLED 
DYNAMICS, SPECIFICALLY DURING THE RETRIEVAL STAGE. MOST OF OUR 
SIMULATIONS HAVE DEALT WITH THIS ASPECT. 
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Figure 1. Geometry of motion. 



SIMULATION PROCEDURE 


EQUATIONS OF MOTION ARE OBTAINED USING HAMILTON'S 
PRINCIPLE, VIBRATIONAL EQUATIONS ARE DISCRETIZED USING RITZ- 
GALERKIN METHOD. THESE DISCRETIZED EQUATIONS ALONG WITH THE 
ROTATIONAL EQUATIONS ARE CONVERTED TO A SET OF FIRST ORDER 
EQUATIONS OF THE FORM x = F(x,e) WHERE e IS A DIMENSIONLESS 
TIME (TRUE ANOMALY) AND x IS THE VECTOR OF STATE VARIABLES, 

THE SYSTEM OF DIFFERENTIAL EQUATIONS IS INTEGRATED USING 
GEAR'S METHOD. 


8 



PHYSICAL DISCRETIZATION: (i) LUMPED MASSES (OR POINT ELEMENTS) 

CONNECTED BY SPRINGS AND DASHPOTS - 
KALAGHAN et al, at SMITHSONIAN 


(ii) FINITE ELEMENTS - KOHLER et al, 

(for ESA) 

MATHEMATICAL DISCRETIZATION: (i) FINITE DIFFERENCE - KOHLER et al. 

(for ESA) 

- KULLA 

(ii) RITZ-GALERKIN METHOD - BUCKENS 

- BANERJEE & KANE 

- BAINUM ET AL, 

- GLAESE & PASTRICK 

- MISRA & MODI 

SEMI-ANALYTICAL PROCEDURES - BERGAMASCHI et al, 
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DEGREES OF FREEDOM 


TETHER ROTATIONS: 


INPLANE (pitch) a 
OUT-OF-PLANE (roll) y 


TRANSVERSE VIBRATIONS 


N 

INPLANE W = L b i (tH i ( y) 

OUT-OF-PLANE u = l J* a(tH (y) 

i=1 


p 

LONGITUDINAL VIBRATIONS : v = l T c t (t)'*' (y) 

1=1 1 1 


THUS THE GENERALIZED CO-ORDINATES ARE 
a, Y j A { , Bj, Cj 

♦ j(y) = 2 SIN (*y/l) 

Tj(y) = (y/l) 1 
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ENVIRONMENTAL FORCES 


GRAVITATIONAL PERTURBATIONS 

DUE TO THE ASPHERICITY OF THE EARTH - J ^ PERTURBATIONS, 
PERTURBATIONS DUE TO THE ATTRACTION OF THE SUN AND THE MOON, 
THESE EFFECTS ARE SMALL, 

SOLAR RADIATION PRESSURE 

CAN BE IMPORTANT AT HIGH ALTITUDES, 

EVEN AT LOWER ALTITUDES, HEATING EFFECT MAY BE SIGNIFICANT, 

ELECTRODYNAMIC FORCE 

FORCE ON A CONDUCTING TETHER MOVING IN THE EARTH'S 
GEOMAGNETIC FIELD. 

AERODYNAMIC DRAG 

THIS IS THE MOST IMPORTANT ENVIRONMENTAL FORCE AT THE LOW 
ALTITUDES (l,E., FOR TSS-J AND TSS-2 MISSIONS). 


11 




WHERE AF IS THE AERODYNAMIC FORCE ON AN ELEMENT HAVING PROJECTED 
AREA Aa 

C D IS THE DRAG COEFFICIENT 

P IS THE DENSITY OF AIR AT THE ALTITUDE OF THE ELEMENT 
UNDER CONSIDERATION 

V VELOCITY OF THE ELEMENT RELATIVE TO THE ATMOSPHERE. 

C D IS ASSUMED CONSTANT IN OUR SIMULATION 


-(H - 


P = P 


h ref^ h o 


-h/h. 


REF 


= P, 


WHERE H = ALTITUDE ABOVE THE SURFACE OF THE EARTH 

OBLATENESS OF THE EARTH IS TAKEN INTO ACCOUNT (AFFECTS h) 
ROTATION OF THE ATMOSPHERE IS TAKEN INTO ACCOUNT (AFFECTS v) 
THE EXPONENTIAL RELATION ABOVE IS AN APPROXIMATE ONE. 

P Q AND H Q ARE DETERMINED FROM THE EXPERIMENTAL DATA SO 
AS TO FIT THE LOWER END OF THE TETHER WELL. 
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PROGRAM 


WRITTEN IN FORTRAN IV TO RUN ON AMDAHL OR IBM mainframe. 
CORE MEMORY REQUIRED - 512K 
SOFTWARE NEEDED - IMSL 


INPUTS - ORBITAL ELEMENTS. 

SEMI-MAJOR AND SEMI-MINOR AXES OF THE EARTH 

AIR DENSITY PARAMETERS AND DRAG COEFFICIENTS 

TETHER PARAMETERS (DENSITY, DIAMETER, YOUNG'S MODULUS, 

DAMPING COEFFICIENT). 

SUBSATELLITE PARAMETERS (MASS, AREA OF CROSS_SECTION) 
INITIAL STATE VECTORS x(O) 

DEPLOYMENT (RETRIEVAL) PARAMETERS 
CONTROL GAINS FOR ROTATIONS AND VIBRATIONS 

OUTPUTS - STATE VECTOR x(e) 

LENGTH, LENGTH RATE 
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SlWULAvTlOM OF 1>F.?LDYKAEMT DYNAMICS 
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Typical deployment dynamics (rotations and inplane orbits 

first modal co-ordinate). 






-V (M/S) L (KM) 

12 . 0 0 2 4 . 0 0 3 6 . 000, 00 25. 0 0 50. 00 75. 00 lpO. 00 



(A) VARIATION OF LENGTH ANO VELOCITY 


FIG ( A-F) DYNAMICAL RESPONSE OURING RETRIEVAL 
FROM 100 KM USING THRUSTERS FOLLOWED 
BY LENGTH CHANGE CONTROL LAW 
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y (DEGREE) CL (DE&R^E) 

- 0.50 0.00 0.50 1.00 - 20.00 0 . 00 20.00 40 . 00 60 . 00 
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-C 2 *10" C, *10' 3 

. 00 0.00 1,00 2.00 3.00 0.00 1 . 00 2.00 3 . 00 4.00 



(C)BEHRVIGUR OF LONGITUDINAL 
GENERALIZED CO-ORDINATES 
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ID) BEHAVIOUR OF OUT-OF -PLANE 

TRANSVERSE HODAl CO-ORDINATES 
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B 2 *10" B, * 1 0' a 

- 0 . 50 0.00 0.50 1.00 - 1.00 - 0.50 - 0.00 0.50 



(E) BEHAVIOUR OF INPLANE TRANSVERSE 
HOOAL CO-ORDINATES 
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0.50 oToO 0.50 1.00 - 5.00 - 2.50 0.00 2. 50 5.00 
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ORBITS 


(F) VPRIflTION OF THRUSTS 
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STRENGTHS AND LIMITATIONS 


STRENGTHS 

“SIMULATION OF GENERAL DYNAMICS INCLUDING TRANSVERSE 
AND LONGITUDINAL VIBRATIONS OF THE TETHER TREATED AS 
A CONTINUUM. 

“CAN HANDLE DYNAMICS DURING CHANGE OF LENGTH QUITE 
RIGOROUSLY. PROVISIONS FOR CONTROL OF RETRIEVAL DYNAMICS. 


LIMITATIONS 

“ROTATIONS OF THE SUBSATELLITE ARE NOT INCLUDED. 
“SOME ENVIRONMENTAL FORCES ARE NOT INCLUDED. 
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FUTURE PLANS 


* EXTENSION TO SPACE STATION APPLICATIONS 

* SIMULATION ON MICROCOMPUTERS 

° SIMULATION OF VIBRATION CONTROL LAWS 
° MULTIPLE TETHER SIMULATION 
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Martin Marietta Simulation 
HIFITS (Hi-Fidelity-Tethered-Satellite) Simulation 

Howard Flanders 

Martin Marietta Denver Aerospace 
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MARTIN MARIETTA HAS SEVERAL PROGRAMS AND SIMULATIONS 
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COORDINATE SYSTEMS AND TRANSFORMATIONS 




COORDINATE SYSTEMS AND TRANSFORMATIONS (CONTINUED) 
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HIFITS SIMULATION - SUMMARY OF EQUILIBRIUM EQUATIONS 



WHERE 
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ABSTRACT DIFFERENTIAL EQUATIONS OF MOTION: 
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INTRO DUCTI ON 

The purpose of this presentation is to discuss the tethered satellite 
system (TSS) portion of the Systems Engineering Simulator (SES). The SES 
is a real time simulator, used at the Johnson Space Center (JSC), to 
evaluate space shuttle operations when there is a man in the control loop. 
The "on-orbit" portion of the SES consists mainly of a mock up of the 
shuttle aft cockpit, complete with functional controls. The operator can 
"fly" the orbiter using the cockpit controls and can respond to simulated 
external scenes and cockpit displays. Other portions of the SES include a 
front cockpit mock up for ascent/descent simulations and a manned 
manoeuvering unit (MMU) for extravehicular activity (EVA) simulations. 

DESCRIPTION 

The on-orbit SES simulates the dynamical interactions between the 
shuttle orbiter, the MMU, and one other payload. (This capability will be 
increased to 4-6 individual payloads in the near future.) All of the TSS 
simulation is controlled from a mock up of the shuttle orbiter aft cockpit. 

This cockpit is complete with hand translational and rotational 
controllers, digital autopilot (DAP) console, radar console, remote 
manipulator system (RMS) console, and the generic CRT 
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consoles/keyboards (MCDS). Status of the tethered satellite can be 
monitored in four ways : 

1 . ) with simulated visual scenes out of the port rear and both 
overhead windows; 

2. ) with simulated camera scenes generated on television 
monitors; 

3. ) with simulated numerical data from the multifunction CRT 
display system (MCDS) display; and 

4. ) with simulated shuttle radar data. 

Manipulation of the satellite control system is done via the MCDS 
keyboard. 

Certain parameters and variables can be output, during a simulation 
run, to both a magnetic tape and an hardcopy printer. The tape may then be 
post processed to generate digital tables or plots. 

During a simulation run, system faults (ie. an RCS jet failure) can be 
generated via an operator’s console. This console will allow faults to be 
input to either the orbiter, the satellite, or the winch system. 
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TETHERED SATELLITE SYSTEM ON SF.S 


62 
















SYSTEM MODELS 


The TSS system is modeled as two free flying vehicles each with six 
degrees of freedom , constrained by an extensible tether with an 
additional degree of freedom. The vehicle translational equations 
include three second order differential equations containing the M50 
spherical coordinates. 

force - mass X translational acceleration (F«mA) - 3 dof 

The vehicle rotational equations include three first order differential 
equations containing the body Euler angle rates and four first order 
differential equations containing inertial quaternions. The Euler angle 
accelerations are integrated to give Euler rates, 

moment « inertia X rotational acceleration (M - la) - 3 dof 

the Euler rates are then transformed into quaternion rates. The quaternion 
rates are integrated to give inertial quaternions which are then 
transformed back to Euler angles for output. 

Shuttle equations of motion : 

^shuttle “ (^ m shuttte^gravity + ^aero + ^RCS + ^tether* ^payload + ^RCS plume) 
a shuttle = I (^gravity grad. + ^aero + ^RCS + ^tether 4 " ^payload 4 " ^RCS plume) 
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Satellite equations of motion : 

^sat * ^gravity + F aero + F APS + F tether + ^shuttle + F RCS plume) 

a sat = ^ ^gravity grad. + ^aero + ^RCS + ^tether + ^shuttle + ^RCS plume) 

The tether, as modeled in the SES, is massless but does include stretch. 
The necessity to keep track of both the stretched and unstretched length 
of the tether introduces one extra degree of freedom. The tether tension 
calculated during the previous time step is input to the natural length 
differential equation. This equation is then integrated twice to get the 
new natural length. The tether stretch (Al) is calculated by subtracting 
the natural tether length from the difference between the positions of the 
orbiter and satellite. The tension is then calculated, by using the stretch 
value, for use on the next time step. In all, there are 13 dof. 

Tether equations of motion : 

- massless, 1/2 of the tether mass is placed at the orbiter CG and 
the other half at the satellite CG 

- the tether mass lumped at the two body CGs is assumed constant 

- the length (A|) acceleration is integrated : 

A| - (R r /J t X- t d + r r f t) 

F t - k Al + 15 km* V A[ 
where : 

R r - radius of tether reel 

J T - moment of inertia of tether reel 

T q - deployer motor torque 
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k - spring constant (AE/I) 

3 - material viscous damping term 
* 

m - generalized system mass 
V A | - time rate of change of stretch 

A - tether material cross sectional area 
E - Young's modulous of the tether material 

Orbiter systems : 

The shuttle orbiter systems are an integral part of SES. These include 
systems such as the digital autopilot, star tracker, navigation, etc. All of 
these systems have been verified by other simulations and experimental 
data. Details of these systems are best found in the shuttle data manuals. 

Rendezvous Radar : 

The SES has a model of the arbiter’s Hughes rendezvous radar. This 
model takes the satellite and orbiter state data and produces output 
similar to what is available on the orbiter : 

1 . ) range / range rate 

2. ) inplane / outplane libration angles and rates 

3. ) line of sight angles and rates 

4. ) etc. 

These data have bias, drift, and random errors introduced to them similar 
to those which occur on the on-board Hughes radar system. 
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At present, the TSS satellite is modeled as a point source with perfect 
reflectance etc. Some preliminary tests show that this may not be a good 
approximation and so some future modifications to the model may be 
required. 

Tethered satellite systems: 

The satellite, on the end of the tether, has three major dynamical 
components as far as the SES is concerned : 

1 . ) mass properties - as previously discussed, the satellite is 
modeled as a free flying rigid body; 

2. ) attitude measurement - the attitude measurements, used for 
control purposes, are assumed to be perfect (ie. no errors are 
introduced); and 

3. ) auxilliary propulsion system - the APS is modeled as perfectly 
aligned cold (N 2 ) reaction jets with step force pulses and perfect 

duration. The libration control jets are controlled by the orbiter 
pilot via the MCDS keyboard while the satellite yaw is controled 
by an onboard control system. The APS fuel usage is monitored 
during usage. 
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ENVIRONMENTAL MODELS 


Gravity model : 

The Earth’s gravitational field may be modeled, in the SES, in two 
different ways, depending on the setting of a digital switch : 

1 . ) spherical Earth (Newtonian gravitation); or 

2. ) oblate spheroid including the first three zonal harmonics (J2, J3, 
J4) and the main Tesseral harmonic (C22, S22) 

Gravity gradient model : 

When an assymetrical body, having one moment of inertia less than the 
other two, is placed in a non-uniform (non-constant) gravitational field, it 
experiences a torque tending to align the axis of least inertia with the 
direction of the gravity gradient. SES is capable of determining gravity 
gradient torques on the orbiter and any free flying payloads. These torques 
may be "turned on" or "turned off" via a digital switch. The TSS satellite 
has essentially symmetrical mass properties and thus, gravity gradient 
torques are negligible. 

Earth atmospheric density model : 

The atmospheric density model used on the SES is the "Bab-Mueller" 
model. This is an approximation to the high fidelity Jacchia atmospheric 
density standard. (See ref. - Jacchia, L.G., "New Static Models of the 
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Thermosphere and Exosphere with Empirical Temperature Profiles", 
Smithsonian Astro-Physical Observatory Special Report 313, May 6, 1970). 

The Bab-Mueller model was chosen because the calculated densities are 
very close to the densities given by the Jacchia model but the calculation 
response is much quicker than for the Jacchia model. The quick response 
time is a fundamental requirement of real time simulations. The 
Bab-Mueller model is described in ref. - (Mueller, A.C., "Atmospheric 
Density Models", Analytical and Computational Mathematics Inc. Technical 
Report ACM-TR-106, June 1977). The classical Bab-Mueller coefficients, 
as given in the referenced report, are run through a Jacchia calibration 
algorithm in the initialization to improve fidelity. Basically the density 
(p) is calculated by : 

P=P 0 e< A + B + c > 
where : 


p 0 - reference density 

A - night time vertical profile exponent 
B - diurnal effect 
C - seasonal (attitude coefficient 

As in the Jacchia model, the sun's orbital position and the effect of the 
solar cycle are considered. In addition daytime/nightime effects are 
considered. 
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Aerodynamic drag : 


The orbiter drag model is a high fidelity model including the orbiter 
attitude effects on drag coefficient and reference area. The drag, pitching 
, and yawing coefficients are obtained by a table look up and linear 
interpolation routine. The roll moment coefficient for the orbiter is very 
small and is considered negligible. 

The tethered satellite is modeled with a constant drag coefficient and 
reference area. This is because the tethered satellite has a very 
symmetric surface profile and the atmospheric drag is independant of 
attitude, (ie. a constant coefficient and reference area). Satellite attitude 
torques caused by drag are not modeled. 

Tether drag is not modeled. 

The atmospheric drag may be "turned on or off" with a digital switch. 
RCS plume impingement : 

During reaction control system operation, the RCS plume may impinge 
upon the surfaces of both the orbiter and the satellite. This can cause 
significant attitude and translational perturbations during proximity 
operations. 
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The plume impingement effects on the orbiter are modeled by curve 
fitting empirical data derived from test firings on an orbiter (ie. curve fit 
constants were developed empirically). These data consider stationary 
orbiter surfaces as well as moving aerosurfaces. 

The satellite plume impingement model can be made up of several 
elements. These elements may be flat, cylindrical, or spherical. If the 
target is more than 5000 ft. from the orbiter CG, no plume impingement 
effects are calculated. The plume flow is modeled as emanating radially 
from a point source, in a direction normal to the nozzle exit plane. The 
plume force/moment is calculated as a function of the jet type, distance 
from the source, magnitude of the thrust pulse, and the impingement 
surface geometry. Details of the jet selection geometry and layout are 
contained in the Shuttle Operational Data Book. Details of the plume 
impingement model are contained in the SES definition document. 

The TSS satellite is modeled as a single spherical plate element. 
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LIMITATIONS 


The TSS simulation has the following limitations : 

- the tether is assumed massless 

- the maximum tether length is 1200 meters 

- the reel motor is not modeled 

ie. Reel torque ■ commanded reel torque 

- external forces do not act on the tether directly 

ie. there is no tether drag 

- a collision between the satellite and the docking ring freezes the 
sim ie. collision dynamics are not modeled 


VALIDATION 

The TSS simulation on the SES was validated by comparing run data 
with data generated by STOCS. STOCS is a non-real time simulation, using 
the TOSS software developed at JSC, with an orbiter digital auto pilot 
(DAP) modeled. 

All of the shuttle systems were validated and are kept valid as an 
ongoing task on the SES. This is done by comparing run data with other 
"off line” simulation data and flight results. 
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FUTURE ENHANCEMENTS 


The following enhancements are being planned : 

1 . ) upgrading the present reel controller to present MMC standard 

2. ) modifying the MCDS CRT display to enhance useability 

3. ) implementing a tether cutting device model 

The following are being considered as possible enhancements : 

1 . ) modeling collision dynamics between the satellite and 
the docking ring 

2. ) modeling the reel dynamics with higher fidelity 
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GTOSS 

GENERALIZED TETHERED OBJECT SIMULATION SYSTEM 


R. PURPOSE fl£ SIMULATION 


The GTOSS software system consists of approximately 410 subroutines representing 
40,000 lines contained in 120 files. Code is constrained to a highly portable subset of 
Fortran 77, with over 700 pages of documentation describing user operation, 
equation derivation, and system software design. GTOSS runs on most computers 
(including: the Macintosh and PC's). GTOSS has been developed under the direction 
of the Avionics Systems Division, Johnson Space Center, NASA. 

GTOSS represents a tether analysis-complex best described by addressing its family 
of modules, designed to be more or less tightly associated as a cooperative whole. 


Tether dynamics: TOSS 


TOSS is a portable software sub-system specifically designed to be introduced into 
the environment of any existing vehicle dynamics simulation to add the capability 
of simulating multiple interacting objects (via multiple tethers). These objects may 
interact with each other as well as with the vehicle into whose environment TOSS 
has been introduced. TOSS is incorporated by adherence to a straightforward set 
of interface rules set forth in the TOSS Interface Control Document. Without small 
motion assumption, and with complete generality, TOSS solves the tether dynamics 
problem relative to a reference point state defined for it by the host simulation. 

Input is designed for easy data identification and entry, as well as to expedite 
parametric studies. Extensive initialization options (such as stabilized gravity 
gradient start-up, Euler angle type selection, etc.) allow user-friendly run setup. 


General tethered system anelysis: GTOSS 

GTOSS is a stand-alone tethered system analysis program, representing an 
example of TOSS having been married to a host simulation. In order to verify the 
TOSS design concept and exercise the TOSS ICD ground rules, it was necessary to 
create a fully representative, yet easily managed simulation into whose 
environment TOSS could be incorporated. The resulting union was called GTOSS 
and has the properties of, and can be viewed as, a system tailored to the 
purpose of examining dynamic behavior of general tethered object 

configurations (space stations, constellations, etc). By contrast, TOSS has also 
been integrated into a Shuttle simulation (to study TSS), with the resulting 
association exhibiting (and rightfully so) the complexity and specificity of an 
Orblter vehicle. The GTOSS host simulation represents a 6 DOF object with 8 
tether attachment points. GTOSS has an executable code size of about 350k Bytes. 
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Input and Initialization for GTOSS (as an entity separate from TOSS) is similar to 
that of TOSS. GTOSS provides output for itself as well as TOSS by invoking 
RTOSS (see below) to generate a Results Data Base to archive solution results for 
display post processing. 


Solution-result archiving: RTOSS 

RTOSS is the Results Data Base (SDB) Sub-system designed to archive TOSS 
simulation results for future display processing. While RTOSS was designed 
primarily to capture TOSS data, it offers a Wild-Card file which can be used by any 
users to capture data of their choice. For instance, GTOSS takes advantage of this 
feature to capture its host simulation data to an RDB for display post processing. 
The modular design of RTOSS requires minimal calls (from the host simulation) to 
invoke the creation and population of an RDB (for instance in GTOSS, this act is 
transparent to the user). At the end of a run, a set of files with unique names will 
have been created containing all pertinent time history data. 

Also inherent in RTOSS are the routines which will extract data from the RDB and 
present it in a form most natural for display post- processing. These extraction 
routines insulate the user from structural knowledge of the RDB, thus rendering 
the user's display software invariant to future changes in RDB design. The RDB can 
be post- processed in a miriad of ways for engineering interpretation, as described 
below. 


Simulation result display: DTOSS 

DTOSS is the first of a growing family of display post processors designed to 
effectively utilize the RDB. DTOSS extracts data from the RDB for extensive 
multi-page printed time history displays. There are currently over 50 different 
display formats to choose from, each of which aggregates data selected to meet 
various display needs. Users are also invited to add new page formats to create 
output for specific needs. 


Simulation result display: CTOSS 

CTOSS is similar to DTOSS, but is designed to create ASCII plot files (as headers and 
data columns separated by delimiters). The same time history data formats 
provided by DTOSS (for printing) are available via CTOSS for plotting. In addition to 
time histories, repetitive tether shape snapshots can be taken in CTOSS. While the 
plot files are optimally configured for existing interactive graphics programs on the 
Macintosh computer, their plot format can be used on other PC's. Since CTOSS 
generates ASCII plot files (which are easily transported between different 
computers), it can be run on any mainframe to generate plot files for a PC or 
Macintosh. 


76 


Simulation result display: ITOSS 

ITOSS differs from DTOSS and CTOSS in two ways. First, it generates output 
display data designed for 3-D animated graphics display of simulation results. 
Second, its output is targetted specifically for the IMI graphics device (a large, high 
resolution, fast display device) . 


Simulation result display: 6ENERRL 


The above RDB post-processors not only represent a significant existing display 
capability for GTOSS, but also function as convenient templates to spawn specialized 
display processors. For instance, hooks are clearly defined, and steps are 
documented to modify the existing DTOSS or CTOSS to add new formats. These 
programs can also serve as a functional boilerplate structure for extracting RDB 
data to be post processed in any fashion you wish. 


B. COORDINATE SYSTEMS flflU DEGREES of FREEDOM 
Coordinate Systems 

The following coordinate systems are used in TOSS: 

A TOSS Inertial Frame The user is allowed to arbitrarily define this frame ; 
however, its relationship to the planet-fixed frame must be explicitly stated in 
the standard TOSS routine invoked to transform inertial frame vector 
components to the planet-fixed frame. This routine (and its inverse- routine) can 
describe an arbitrarily complex relationship between inertial and planet-fixed 
frames. As delivered, the TOSS inertial frame is defined as one which is aligned 
with the planet-fixed frame at zero simulation time. 

B. Planet-fixed Frame. This is typically an earth-fixed frame, and also the one in 
which all planetary environment calculations are defined. The user is allowed to 
arbitrarily define this frame, however, its relationship to the inertial frame 
must be explicitly stated in the TOSS routine which transforms planet frame 
vector components to the inertial frame, and environment calculations must be 
consistent with its definition. As delivered, GTOSS environment routines assume 
an earth- fixed frame, the +X axis of which is presumed to pass through the 
Greenwich meridian, +Z axis through the geocentric North pole. 

C. Topocentric Frame. This is a frame aligned along local spherical longitude, 
latitude, and radius vector to the planet center It is used for state initialization 
options and result interpretation. 

D. Orbital Frame. This frame is defined by the kinematic state of a point, and is 
similar to the topocentric frame except the local longitude/latitude vectors are 
aligned to the current plane of Keplerian motion. Any TOSS entity can have an 
associated orbital frame, which is used for state initialization options and result 
interpretation. 
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E. Body Axis Frame. A body-fixed frame (and a body station reference point) is 
associated with each 6 DOF object. The frame and body station are arbitrary, 
but must be consistently used for defining all body attributes (CG location; 
tether attach-points; aerodynamic reference point, etc). This frame is the 
reference for body attitude interpretation. 

F. Tether Frame. Each finite tether has its own tether frame. The X axis of this 
frame is aligned along the line of sight between the tether's attach points. The Z 
axis is orthogonal to the first, but lies in the orbital plane (of a preferred 
kinematic point). The Y axis is defined to complete the triad. This frame hosts 
the tether dynamics coordinate solution, and is used for both initialization 
options and result interpretation. 


Degrees of Freedom (DOF) 

Many of the DOF are modifiable via procedures included in the manuals The 
nominally delivered GTOSS configuration provides: 

1, Up to 9 bodies (each with up to 6 rigid- body DOF's). 

2. A 3 DOF particle dynamics option to eliminate rotational dynamics overhead 
(for efficient study of overall topological behavior). 

3 Concurrent simulation of up to 25 different tethers. 

4. Tether inter-connection in any conceivable fashion at up to 8 attach-points 
per body. Each attach point is specified by its coordinates in an arbitrary 
body axis frame 

5 No practical constraints on tether/attach point connectivity. 

6 Tether dynamics simulation as your mixed choice of 

• Massless models (ie. linear springs and dash pots) 

• Modal Synthesis finite models (on a tether-by-tether, and axis-by-axis 
basis, up to 15 modal coordinates per each of 3 axes, with the 
evaluative resolution of each generalized force type specifiable at up to a 
maximum of 30 uniform spatial intervals). Tension is optionally 
evaluated either in terms of material strain, or Lagrangian multipliers 

• Point Synthesis (ie. Bead) type finite models (on a tether- by- tether 
basis, up to 48 collocation points per tether; with 3 DOF per point). 
Tether can break on tension limits, or be severed at multiple points. 

7. Five length-rate, and five tension-profile, data-driven deployment scenarios: 

8. Five power-profile (amp-limited} data-driven, power generation scenarios: 

• All deployment and power scenarios can be arbitrarily assigned to any 
one or more tethers (a tether also has its own scenario scaling factors). 
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C. ENUIBONMENTRL MODELS 

Planetary Enulronmant 


Subroutine calls to the planetary environment models are standardized within TOSS 
so that any level of environmental sophistication (you have available) can be easily 
incorporated into TOSS. Calling arguments to an environment model are: a 
sophistication-level flag; and, a position state in the planet-fixed (earth) frame. 
Vector results are returned in the planet-fixed frame. Currently delivered with 
GTOSS are the following: 

o Gravitational Field Model: Earth; Inverse square central force field, also an 
Oblateness model with 2 anomaly terms. 

o Globe Shape Geometry Model: Earth; Spherical globe. 

o Atmospheric State Model: Earth; 1976 Standard earth density, speed of 
sound, and temperature model (Z% accuracy to 1000 KM). 

o Atmospheric Kinematics Model: Earth; Rotating atmosphere (with local wind 
perturbations currently zero). 

o Magnetic Field Model: Earth; Tilted, shifted, vector dipole. 

o Inertial Frame Model: Planet centered inertial point, with inertial frame 
aligned to planet-fixed frame at zero time. 


Entity Attribute Environment 


As delivered, all TOSS objects can experience simple aerodynamic drag. Tethers 
experience distributed aerodynamic lift and drag. Tethers carrying electric current 
experience electromagnetic forces. 

TOSS is designed to facilitate incorporation of any degree of entity attribute 
simulation. 


0. SYSTEM MODELS 

The design philosophy of TOSS has been to provide a useful array of built-in system 
simulation features such as data driven scenarios for: control forces and moments, 
tether deployment, and electromagnetic power generation; as well as default 
aerodynamics options; etc. In addition, TOSS supports arbitrarily complex 
simulation of mass properties, control systems, aerodynamics, tether deployment, 
and power generation, etc. through its documented user-interface data structures, 
modularized code, and logical hooks which invite and assist in user modification. 

Of course, when TOSS is incorporated into a host simulation, all the system models 
present in the host then function in the tether environment provided by TOSS. 
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E. L1MITHTI0NS 


The following summarizes the currently known limitations to the use of 
GTOSS/TOSS (also see section below entitled: What should be done next). 

a. Tether Frame Definition: The orientation of the Tether Frame is undefined if 
attach point line-of-sight becomes exactly perpendicular to an associated orbital 
plane. Near this state, frame orientation rates can become large, inducing 
large apparent rates of change of finite tether coordinates. This can be avoided 
in the Point synthesis model by integrating coordinates in the inertial frame 
(thus using the tether frame for interpretation only). To date, this has not 
been a practical restriction (due to the nature of engineering applications); If 
required, this restriction can be removed. 

b. Modal Synthesis Finite Tether Model: This model is only valid during motions 
for which the distance between attach points is greater than the deployed 
tether length (a state of intrinsic tension). Furthermore, cases in which 
cyclic slack/taut states are a significant element of behavior may not be 
simulated well. In short, while the Modal synthesis model has some 
advantages, the Point synthesis (bead) model is significantly more robust and 
should be used to establish truth datum. 

c. Stabilized Gravity Gradient Initialization: Currently this feature applies only to 
simple chains, limited to 3 objects and 2 tethers (as a mixture of massless and 
Point synthesis models). Multiple disconnected chains are allowed. 

d. Certain Environmental Effects (for example, solar pressure effects) are not 
incorporated into TOSS. 


F. URL I POTION MDMflflS 

Validation of TOSS/GTOSS involves three distinct areas: Classical techniques (for 
those cases which are simple enough to be described by classical closed form 
solution); Comparative techniques (for those cases which defy classical 
substantiation); and Official techniques (to verify site installation and validate 
evolutionary changes to GTOSS). 


Classical Uertflcatlon 

Offical solution verification of TOSS has been accomplished from many directions. 
First, the simulation of the rigid body TOSS objects are verified via classically 
known solutions to Euler's rotational equations and Newton's law. Overall 
translational dynamics is further verified against classical Keplerian motion. 
Massless tether dynamics are verified against known gravity gradient as well as 
bolo type motions. Finite tether wave propagation and shape verification has also 
been accomplished against classical string theory. 
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Comparative Uerification 


Finite tether dynamics have been verified by comparing both the TOSS Modal 
Synthesis (MS) and Point Synthesis (PS, ie. bead model) against each other and 
against a bead model independently developed by NASA JSC. These comparisons have 
addressed aerodynamic response, wave propagation, and slack/taut behavior. 
Electrodynamic finite tether response has been officially compared only internally 
to TOSS (between the MS and PS models). 

There are at least 5 known, actively used, installations of GTOSS. When GTOSS is 
installed, users invariably compare its results to independent tether solutions 
previously verified by the installation. To date, there are no unexplained anomalies 
in result comparisons of this type. 


Installation Delineation 


Delivered with GTOSS are 18 different input run decks (along with corresponding 
official results) which test all aspects of GTOSS. A user site compares the output of 
the newly installed GTOSS to these official results to verify site installation. 


6. IDHHT SHOULD BE DONE NEXT 


The following features are either planned, or recommended for future GTOSS/TOSS 
development: 

1. Non-uniform point spacing for the Point Synthesis (PS) finite tether model: 
This will allow more efficienct use of degrees of freedom for solving certain 
type finite tether problems. 

2. Non-uniform material properties for the PS finite tether model: This will allow 
simulation of tethers which are purposely constructed of different materials as 
a function of length. 

3. Auto- transition (on user defined criteria) from finite to massless to finite 
tether models: This will allow full mission simulation continuity without the 
large computation overhead associated with very short tethers. 

4. Improved aero force model for finite tethers: The current model does not 
represent the consensus standard for distributed air loads on tethers 

5. Flexible boom attachment simulation. This will allow a certain amount of 
attach structure flexibility without full involvement in simulating body 
flexibility of a TOSS object. 

6. While the latest version of TOSS reflects an emission-end-sensitive plunging 
aspect of tether deployment (under small strain), the author does not admit to 
knowledge of an ultimate expression for deployment kinematics or dynamics. 
TOSS will be continually improved in this area as understanding is gained. 
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H. £M £ STUDY EHRMPLES 
Examples of Host/TOSS Integration 


There are two different instances of TOSS having been successfully introduced into 
a host simulation environment: the first of these is GTOSS (described above), the 
second is called STOCS (Shuttle Tethered Object Control Simulation). The host 
simulation for STOCS is a full engineering fidelity flight control simulation of the 
Orbiter. STOCS is being used to perform mission verification tasks for the TSS 
experiment (under the responsibility of NASA JSC). STOCS is an excellent example 
of complex control systems (TSS on-board control) becoming associated with a TOSS 
object (which represents the TSS), as well as the introduction of specialized tether 
deployment systems (TSS deployer) into TOSS. 


Examples of TOSS/GTOSS Application 

TOSS has been used to study enumerable engineering applications of tethers. Some 
of these are: 

1. TSS mission study (STOCS/TOSS, 2 rigid-bodies, both massless and bead model 
tethers). 

2. A spinning, orbital station ( 6 rigid-bodies, 15 massless tethers). 

3. Gravity gradient/Gyroscopic orienting spinning dumbell (2 particles, 1 massless 
tether). 

4. Electrodynamic day/night power gcneration/orbital boost with tethered counter 
balance (3 particles, 1 bead model and 1 massless tether), 

5. Real-time Shuttle Engineering Simulator (SES) verification (comparison runs 
made with both GTOSS and STOCS). 

6. Planetary exploration maneuvers using a slingshot mechanism. 

7. Electrodynamic pulse maneuvers (2 particles, 1 bead model tether). 

8. Space Station docking devices (3 particles, 2 massless tethers). 

9. Orbiting, spinning carousel (3 rigid-bodies, 1 massless, 1 modal synthesis, 
and 1 bead model tether). 

10. Gravity gradient stabilized orbiting platform (3 rigid-bodies, 1 modal synthesis 
and 4 massless tethers). 

11. Simulated bead model by chain-connecting 9 TOSS rigid-bodies with 8 massless 
tethers (for bead model verification study;. 

12. Comparative solution verification studies between GTOSS and an independent 
bead model simulation (wave propagation, aerodynamic response, symmetrical 
slack/taut gravity gradient behavior of a system exhibiting TSS physical 
properties) 
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| GTOSS TERMINOLOGY || 

TOSS II? Tethered Object Sub-System 

GTOSS ItF 3 Generalized Tethered Object 

Simulation System 

nmm KS 3 Result Data Base sub-system 

nmm Display print ROB post-processor 
CTFGD§)S§ CF 3 Cbart/Graphics ROB post processor 
imm kf* i mi graphics ROB post processor 
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INTERFACE CONCEPT 





SVSTEM 

INTERFACE 



INTERFACE DATA STRUCTURE 

• REF. PT. STATE UJR/T INER. PT. 

• ATTACH POINT KINEMATIC STATES 

• ATTACH POINT FORCES/COUPLES 


cf- No Orbital State Assumptions 
No Small Motion Assumptions 
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GTOSS OPERATION 




GTOSS Input 
file INDATA 


means 
HSCII file 



Optional CRT 
interactive input 

Optional CRT 
run Monitoring 


Quicklook output 
file OUTDATA 


Results Data Base 
file set run / 

i 

Results Data Base 
file set run 2 


Results Dat8 Base 
file set run N 


HRCHIUED ROB’S 


DTOSS Input 
file INDIS 


DTOSS 


Multi -Page Print 
Facility output 
file OUTDIS 


Optional CRT 
interactive input 

Optional CRT 
run Monitoring 


(Similar RDB 
processing 
for CTOSS 
and ITOSS ) 
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9 Bodies (3 or 6 DOF) 


25 Different Tethers 


8 Tether attachment points per body 


Totally General Connectivity 


Tether models, your mined choice of: 



Massless 

Modal Synthesis (3-D) 


Up to 15 Modal Coordinates (Legendre) per oh/s 
Up to 30 Generalized Force eualuation interuals 
Tension uia strain or Lagrongion mu/tipiers 



Point Synthesis (3-D) 


Up to 50 collocation points per tether 

3 DOF per collocation point 

Tether seuer (on time) or break (on tension) 


Data-driuen, arbitrarily assignable scenarios 

5 length-rate deployment scenarios 
5 tension control deployment scenarios 
5 pouter profile (amp limited) generation scenarios 
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Input to all TOSS enuironment routines: 

1. Fidelity leuel flag 

2. Simulation Time 

3. Position in Planet-fined frame 


Inertial Fram^Model | EFTEI, EITeT 

To eualuate this enuironment, use standard routine 


Grauity Model | GRRU 
GlobeShap^Mode^j EE0D~ 
Magneti^neldModel | GRUSS 
fltmospherl^tate Model j) nTM0S 


Atmospheric Kinematics Model | WINDS 
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Tether-Frame is Un-defined uihen 
line of attach points is 90 degrees 
out of the orbital plane 


$ Modal Synthesis model ualid only 
for a state of intrinsic tension 


*=> Modal Synthesis model dubious 
for cyclic slack- foot states 

9 Stabilized Grauity Gradient start-up 
ualid only for simple chains of up to 
3 bodies and 2 tethers (mixed bead 
and/or massless) 
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MMrVHrVHrVH J lMr4MrVHrVHrVHrV^ 


CylCLRSSlCRL UERIFICRTIDN 


♦ Keplarian/Neiotonian behauiour 


• Classical String Theory 


Cj^|COMPRRRTIUE UERIFICRTION 


Participants: 

• JSC deueloped bead model 

• GTOSS bead model 

• GTOSS modal synthesis model 


Study parameters: 

■ TSS Type Parameters 

• Hero Response 

■ Transuerse Ulaue Response 

• Symmetrical, Un-forced Response 


EP’IOFFICIRL UERIFICRTION 


K^IUN-OFFICIRL UERIFICRTION 
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TOSS/GTOSS RPPLICRTION 

• TSS/Orbiter Compatibility (2RB, 1 TH) 

• Spinning, Orbital station (6RB, 15TH) 

• Grauity gradient/Gyro dumbell (2 P, 1 TH) 

• Day/Night electro pivr gen (3 P, 2 TH) 

• SES uerification (GTQSS/STQCS) 

• Planetary exploration studies 

• Electro pulse maneuvers (2 P, 1 th) 

• Space Station docking deuice (3 P, 2 TH) 

• Orbiting spinning carousel (3 RB, 3 TH) 

• Gravity gradient platform (3RB, 5 TH) 

• 9 Body equivalent bead model (8 TH) 

• Uerification studies (20 - 50 beads) 
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WHIRLING DERVISH 



91 




| FUTURE GTOSS PLANS [ 

C^> Non-uniform Point Spacing 

C^> Non-uniform Tether Properties 
t£> fluto-phase transition 
Cj>lmproued aero-force model 
C^> Flexible boom simulation 

l^>lmproued deployment fidelity 
C^> Expert, Friendly User interface 
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(FT) 


BEAD MODEL V8 MODAL SYN AERO RESPONSE 



OJO 0.2 0.4 0* 0J8 1 JO 


NORMALEED TETHER LENGTH (300k FT) 


h 


( 

a 

2 

ill 


BEAD MODEL WAVE PROPAGATION (48 Beads) 
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(snon) JU3 tauai miuM/winuMi 


DRY/NITE 


6ENERHTI0M (65k LENGTH) 






Tethered Satellites Simulation 

John R. Glaese 
Control Dynamics Company 
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CONTROL DYNAMICS COMPANY 
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Dr. John R. Glaese 
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X vs. Y plots 


Simulation Coordinate Systems 
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(Local Vertical) 


Partial Differential Equation (PDE) 
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Dynamics Equations of Motion 
Boundary Conditions 
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DISCRETIZATION 



104 


DISCRETIZATION 

(Continued) 
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TENSION CONTROLLED DEPLOYMENT 
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Amiga Graphics Worksatation Development 

Work Sponsored by NASA/MSFC 

Objective: Develop tether simulation graphical 
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150 kg satellite 
upward deployment 
24 km deployment 
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TENSION (N) vs. TIME (sec 
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Skyhook Program 
David A. Arnold 

Smithsonian Astrophysical Observatory 
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SKYHOOK PROGRAM 


A. PURPOSE OF SIMULATION 

• GENERAL PURPOSE PROGRAM FOR SIMULATING THE DYNAMICS 

OF TETHERED SATELLITE SYSTEMS 

• FEASIBILITY AND DESIGN STUDIES 

• STUDY SCIENTIFIC AND OTHER APPLICATIONS OF TETHERED SYSTEMS 


hkchhng page BLANK FHJ<ED 
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B. COORDINATE SYSTEMS AND DEGREES OF FREEDOM 


• RECTANGULAR GEOCENTRIC COORDINATE SYSTEM 
. POSITION AND VELOCITY OF EACH MASS 

• TEMPERATURE OF THE WIRE 

• ORIENTATION OF THE SUBSATELLITE (DIRECTION COSINES) 

• ELECTRIC CHARGE ON EACH MASS 
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C . ENVIRONMENTAL MODELS 


• SPHERICAL HARMONIC REPRESENTATION OF EARTH'S GRAVITY 

• LUNAR & SOLAR GRAVITY PERTURBATIONS 

• ATMOSPHERIC DRAG FORCE 

• SOLAR, EARTH INFRARED, AND ATMOSPHERIC DRAG HEATING 

• EARTH'S MAGNETIC & ELECTRIC FIELDS 
. IONOSPHERIC DENSITY 


D . SYSTEM MODELS 


• BEAD MODEL OF THE TETHER WITH VISCO-ELASTIC FORCES 

• AREA TO MASS RATIO OE SUBSATELLITE, SHUTTLE AND WIRE 
. TORQUES ON SUBSATELLITE 

• THERMAL ABSORPTIVITY, EMISSIVITY, AND HEAT CAPACITY 

OE THE WIRE 

• THERMAL EXPANSION COEFFICIENT 

• CAPACITANCE AND RESISTANCE OF THE WIRE 

• CHARGE COLLECTION MODELS OF THE WIRE, SHUTTLE, AND 

SUBSATELLITE 

• ELECTRODYNAMIC FORCE ON THE WIRE 

• THRUSTER MODELS 

• VARIOUS REEL CONTROL ALGORITHMS 
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E. LIMITATIONS 


• COMPUTER TIME AND MEMORY REQUIREMENTS INCREASE RAPIDLY 

WITH THE NUMBER OE MASS POINTS 

• RESOLUTION OF WIRE DYNAMICS IS LIMITED BY NUMBER OF 

MASS POINTS 

• RESOLUTION VARIES DURING DEPLOYMENT & RETRIEVAL AS 

BEADS ARE ADDED OR SUBTRACTION 

• INDUCTIVE EFFECTS NOT INCLUDED IN ELECTRODYNAMIC MODEL 
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F. VALIDATION METHODS 


• COMPARISON OF SIMULATIONS WITH ANALYTIC SOLUTION FOR 
THE BEHAVIOR 

. COMPARISON OF RESULTS WITH DIFFERENT NUMBERS OF MASSES 
. COMPARISON OF RESULTS FROM DIFFERENT SIMULATION PROGRAMS 
. EXAMINATION OF INTERMEDIATE CALCULATIONS DURING DEBUGGING 
AND COMPARISON WITH HAND CALCULATIONS 
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G. WHAT SHOULD BE DONE NEXT 


. MODIFY SKYHOOK TO MAINTAIN RESOLUTION AS WIRE LENGTH 
CHANGES 

• COMPARE RESULTS OF BEAD MODEL OF WIRE WITH MODAL 
REPRESENTATION 
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Figure 6f 
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SLACK3 


Gordon E. Gullahorn 
Smithsonian As trophysical Observatory 
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SLACK 3 
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THE MODEL 
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LINEARIZED GRAVITY GRADIENT AND CORIOLIS FORCE 
ATMOSPHERIC DRAG 
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SLACK3 INPUT PROMPTS 

These prompts define the parameters available for 
specification by the program user. 

THE FOLLOWING PARAMETERS AND FLAGS HAVE DEFAULT VALUES AND 
MAY BE CHANGED BY ENTERING THE DESCRIPTOR IN THE FIRST COLUMN. 


ORIENTATION 
TETHER. PROP 
TETHER. ANGLE 
TETHER. TENS 

DISCRETIZE 

ORBIT 

DEBUG 

SEG . LENGTH . RANDOM 

SEG. SLACK 

SEG . DIR 

DRAG . FUDGE 

BOOM. PROP 

BOOM. ANGLE 

BOOM. PHASE 

ROT. AXIS 

ROT . RATE 

OFFSET 

REL.ABS 

BURNS 

CUT. FREE 

RUN. TYPE 


UP/DOWN TOGGLE IS: UP 

1 DENSITY 1.40 E 0.14D+12 CV 0.30D+09 DIAM 0.30 
PITCH 0.00 ROLL 0.00 
ORIG LENGTH (KM) 20.0 SS MASS (TONS) 0.50 

TENSION FUDGE 1.00 
LENGTH RATIO, ATTACHMENT/TIP: 0.50 
QUADRATICNESS (0-LINEAR TO 1): 0.50 
HEIGHT (KM): 295.0 ATM. DENS.(G/CC): 0.27D-13 
DEBUG OUTPUT TOGGLE: F 

% SEGMENT LENGTH RANDOMIZATION: 5.0 

SLACK FACTOR: 0.95 % RANDOM: 5.0 

SEGMENT DIRECTION RANDOM : 5.0 

BALLOON DRAG FUDGE FACTOR: 0.75 

El, LENGTH, WEIGHT (NASA UNITS): 0.13D+09 849.0 106.2 
PITCH, ROLL: 0.0 0.0 

PHASE AT T*0 : 0.0 

DEFINING ANGLES, PITCH & ROLL 90.0 0.0 

RATE & ACCELERATION (DEG, SEC): 0.00 0.00 

X, Y, Z (METERS): 0.00 0.00 0.00 

RELATIVE/ABSOLUTE OUTPUT TOGGLE: REL 
0 THRUSTER BURNS. (MAXIMUM 5 ) 

CUT TETHER FREE AT T- 9999.9 
BREAK 


THE FOLLOWING INPUTS ARE ALWAYS REQUIRED: 

ENTER RANDOM SEED, ODD INTEGER < 2147483647 
NUMBER OF TETHER SEGMENTS? 3 TO 50 

LENGTH AT TIME OF REEL JAM (KM)? 

(The program now computes: ) 

A CRUDE TIME SCALE FOR RECOIL, NEGLECTING GRAVITY 
GRADIENT, IS: 1329.333973151609 

HOW LONG (TETHER-SECONDS) TO COMPUTE? 

HOW OFTEN (SECONDS) TO OUTPUT? 

SET A LIMIT ON # OF BOUNCES (INTEGER): 
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Figure 2.4.4. A* for Fig. 2.1.1, except: Thrusters directed along-orbit 

cutoff after 20 second burn, total run of 1200 seconds. 


38 





Figure 2.4.3. As for Fig. 2.1.1, except: Thrusters directed perpendicular to 

orbital plane, cutoff after 20 second burn, total run of 1075 seconds. 
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TIME (SFO 






ORDINAL PAGE IS 

Of POOR QUALITY 


SIMULATION A^_ SIDE-VIEW 






'iW M 'NM. M M VlM H -JMM M >MM) VlM M H B >«M M > «M « >MC M -W«.» -^IM M Vn».« 

- Expended Scale, Side-View Plotted Every AO Sec for 7200 Sec. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 

SIMULATION JL. SIDE-VIEW 




-Me M M IVa oe ~*oi 
IN-FIRNE C0«FONEMT mETEASl 


W.M owl I 



~ow.«“ow.« oura -UM.M -IM.I ->«•" Amt.m 

- Expanded Scale, Side-View Plotted Every 40 Sec for 7200 Sec. 
Longitudinal + In-Plane Librational/Lateral Dampers On. 
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SIMULATION A AND JL. FRONT- VIEW 


MASS SYMBOL 

1 O 

2 a 

3 + 

4 X 



t iU.m iW.m iS.m 3T«i «Tm Tm "T5T5 -W'm mn.n *5w!n -Wm ThTn H5 ' Xui.m ''4m7m 3wT 



m.m 4m. m 4*Tm im ao -''m oo ■'»« * 
OUT-Of-PlflNE COMPONENT INETEAS) 



■TM.M -IM.N -'iMa.M -I1M.M -IIM.N -IJW.M -tlM.H -I1N.M -TlM.M •'tMB.M -'iMt.M -IU«.M -TiM.M 


- Expanded Scale, Front-View Plotted Every 40 Sec for 7200 Sec. 
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PLAN OF FUTURE ACTIVITY 
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PLATFORMS, FURTHER DEVELOPMENT OF CONTROL/NOISE -ATTENUATION 
STRATEGIES, INVESTIGATION OF FORCED DYNAMICS OF 4-MASS 
TETHERED SYSTEM. 



Elasticity Effects in Tether Dynamics 

Silvio Bergamaschi 
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OBJECTIVES OF THE ITALIAN PROPOSAL 


The overall purpose is to understand the relevance of tether elasticity 
in TSS dynamics, because: 

long tethers have never been flown; experience is lacking 

high frequency dynamical noise might preclude future 
experiments ( gravity gradiomet ry , atmospherics ) 

- non-norainal situations can occur where elastic effects are of 
primary importance ( tether slackness, impact with a 
micrometeoroid) 


Thus, the work plan is: 


To investigate natural frequencies of elastic vibrations of TSS 

To study amplitude of response to excitations 

To analyze data from gyros and accelerometers mounted on the 
satellite 

To check theoretical and experimental results 
If needed, to improve the simulation models 


THE MATHEMATICAL MODEL 


Assumptions : 

Space Shuttle in circular, unperturbed orbit. Earth oblate* 
ness, atmospheric and electrodynamic drag neglected 

Perfectly elastic tether. No material damping included 

* Point mass satellite 

Constant tether length 

Constant tether diameter 

Ambient temperature changes ignored 


Method 

Lagrangian formulation adopted 


°HECEDiNO P^GE 


blank not filmed 
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TRIAL AND EROR PROCEDURE TO COMPUTE THE EIGENVALUES 


To assume a tentat ive value for q 

To integrate numerically the resulting c equation (in 
nondiraensional variables) in the interval 0-1 with c(o)=0 and 
c"(o)=l as initial conditions 

To check if the boundary condition at the satellite end of the 
tether is satisfied. If not, to reflect the sequence with an 
increased |q| value. 
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FIRST THREE MODES OF OUT OF PLANE MOTION 


TETHER STRUCTURAL DAMP I N G 


Purpose: To investigate transients in longitudinal vibrations 

Assumptions : 

- coupling with orbital motion and lateral vibrations 
neglected 

experimental data from Martin Marietta used 


Dynamical equation: 





oN 

+ C 

oV> 


(oufre 


Boundary conditions: 


eft y (.«,» = -/*• j 


Initial conditions: 


°) s K* j (*. o) - o 
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SYSTEM DYNAMICS • RIGID TETHER 


Purpose of the investigation : 

To analyze system stability during one elevator transfer from 
the Space Station to the SATP. 

To evaluate orbital perturbations on the Space Station, 

To devise velocity control laws, if needed. 

Assumptions : 

Tether elasticity neglected 

Characteristic dimensions of Space Station, elevator and SATP 
negligible with respect to tether length 
In plane motion only (Coriolis force) 

- No external perturbations taken into account (Earth oblateness, 
atmospheric drag, etc,) 


Method : 


Lagrangian formulation of the equations of motion 


THE SYSTEM 


Retrievable platforms tethered to the Space Station in order to: 

- conduct experiments in a dynamic noise force environment (not 
attainable on the Space Station) 

Elevator envisaged as: 

platform for microgravity experiments at system center of 
gravity 

servicing module for the SATP 

Tether intended to: 

sustain SATP gravity gradient force 
act as structural damper 

transfer power and information (umbilical concept) 
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Free notion 
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FREE MOTION 
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CONTROLLED TRANSFER 


Control laws affected by: 

Coriolis force, with impact on average translation velocity 

- Transients at beginning and end of each maneuver, because of 
singularities of the dynamical equations at both ends of the 
tether 

linear Control : 

In some cases numerical instabilities at the very beginning of 
maneuvers 

Exponential Control : 

At the beginning: 

- At the end: 

- Smoother motion 

- Reduced time for a transfer 


= v. WU*) 

d ( - i/o (4 > 4 *) 


Slightly increased perturbations on Space Station orbit 
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_J _ 2 PERTURBATIONS 


Purpose: To evaluate tether motion around the local vertical, forced by 

Earth oblateness 

Assumptions: 

tether elasticity neglected 

first order theory used for the changes of the oscillating 
elements 

- three-dimensional motion simulated 
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TSS SUBSATELLITE ATTITUDE DYNAMICS 


AND CONTROL LAWS VERIFICATION PROGRAMS 
F. Venditti 

Aeritalia - Gruppo Sistemi Spaziali 
Torino, Italy 


Abstract- This paper deals with the presen- 
tation of a dynamic model of the Tethered 
Satellite System and of the relevant 
simulation program, developed in order to 
provide the dynamic analysis support for the 
design verification of the Subsatellite 
attitude control, A special care has been 
spent in the Satellite attitude dynamic 
analysis and the model has been specifically 
conceived to this aim. 

The way in which the simulation results can 
be utilized for the verification and testing 
of the attitude control is also presented. 

1 - INTRODUCTION 

The design of the Attitude Measurement and Control 
Subsystem (AMCS) of the TSS-Subsatellite for the electrody- 
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namic mission has been completed and now it is going to be 
verified and tested. A representative dynamics of the system 
is one of the main items that has to be implemented for the 
design verification phase and in such an optic a dedicated 
program has been developed. 

The aims of this program are: 

- to investigate the dynamics of the system, with a 
particular interest for the Satellite attitude behaviour, 
and to identify potentially critical areas for the AMCS 
performance verification; 

- to provide realistic input data for the verification 
and testing phases. 

To this purposes a model representing the overall 
Orbiter-tether-Sateliite system has been considered and 
developed, while the adoption of a simplified, provided 
representative, AMCS model is enough for the goal. 

On the other side, to the AMCS performance verification 
aim, a complete and accurate model is needed for the AMCS 
itself, requiring a small integration step (128 msec), tied 
to the high on-board data handling sampling rate. The need 
for the implementation of a simplified dynamic model (3 
satellite rotational d.o.f + the possible elastic ones) 
arises then, in order to avoid unaffordable CPU times. 

The "dynamic-simplified" model can be called an 'open 
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loop' approach because it assumes that the overall system 
dynamics impacts on the Satellite attitude, while the fact 
that the attitude is controlled has a negligible effect on 
the overall behaviour, this resulting in the tether tension 
at the swivel point which does not depend on the Satellite 
attitude. This seems to be a reasonable hypothesis, to be 
anyhow validated by means of the "AMCS-simplif ied" model, 
which, in turn, can be called the 'closed loop' approach, as 
it accounts for the coupling between the satellite attitude 
and the overall dynamics, while being not burdened by low 
integration step requirements. 

The simulation results of the 'closed loop' model will 
then be used as input data for the "dynamic-simplified" 
model runs; in fact, in order to provide an adequate capa- 
bility to represent the system, the "dynamic-simplified" 
model requires, for each run, the force time history at the 
tether attachment point and the firing sequence of the 
Satellite thrusters not driven by the AMCS (in plane and out 
of plane thrusters) in input. 

The results of the 'closed loop' simulations will also 
be used for the integration and testing phase: the AMCS 
hardware and software will work, by means of gyros and 
sensors stimulators, in closed loop with a simulated real- 
time-computed system dynamics. The same 'open loop' model 
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as for the verification analyses will be used to provide the 
necessary input to the stimulators. 

At the moment, the 'closed loop' model has been 
implemented and the program is undergoing integration tests, 
while the first simulation runs will be performed in a short 
time. 

Hereafter the 'closed loop' model and program are 
described. 


2 - MODEL DESCRIPTION 
2.1 - Configuration 

The electromagnetic mission will be supported by the 
Shuttle in an approx 300 km orbit; the Satellite will be 
deployed upward up to 20 km from the Orbiter; the Satellite 
mass will be approx 500 kg. 

Only the Satellite yaw axis is controlled, stabiliza- 
tion on the other two axes being assured by the gravity 
gradient. During on-station phase, the Satellite spins about 
the yaw axis and the spin rate is controlled (1 rpm) . During 
retrieval, the yaw angle is controlled in such a way to 
allow for the in plane and out of plane thrusters to fire in 
the correct direction for libration control. During the on 
station phase, two deployable booms, part of the scientific 
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experiment, will be extended to approx 4m. 


2.2 - Tether Model 

Some preliminary analyses showed different kinds of 
behaviour of the system depending on the tether length. This 
leads to different representations of the system, depending 
on the deployed length. In figure 1 the main system frequen- 
cies are given as a function of tether length (£). 

The frequencies of the Satellite oscillations about 
its roll and pitch axes decrease as H decreases; the gravity 
gradient stabilizing effect for these axes also decreases, 
support in tensioning the tether being provided at short 
distances by in-line thrusters. 


Fig. 1 
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The longitudinal frequencies of the overall system due 
to tether elasticity decrease as i increases. It can be seen 
that the axial frequencies may couple with the pitch/roll 
oscillations. The number of modes to be taken into account 
depends on the tether length: at short distance only the 
fir3t axial mode is significant, while the others could be 
neglected, leading to a reduction of the integration step 
without affecting the system behaviour representation; at 
the longest tether extension also the second mode has to be 
considered. 

The string vibration too may have a significant 
coupling with the Satellite oscillation at short tether 
length. 

These considerations lead to the need of considering 
in any case the tether mass and elasticity: a lumped mass 
approach for describing the tether was chosen, where the 
number of masses is selectable as a function of tether 
length. 

The tether can be modelled as a massless spring 
connecting the Orbiter and the Satellite; in such a way only 
the first axial mode is simulated. To include the other 
modes, the tether is described by mass points connected by 
springs whose stiffness depends on the undeformed tether 
length between adjacent bodies. 


194 


During deployment and retrieval the tether mass chan- 
ges: this is represented by changing the mass of the mass 
point closer to the Orbiter; the other mass points maintain 
their mass. The nominal (unstressed) tether length between 
the Orbiter and the closer mass point also changes (and so 
its stiffness does), according to tether reeling, while the 
unstressed tether length between the other masses remains 
unchanged . 

The number of mass point cannot be changed during the 
simulation in the current version of the program. The fixing 
of this number at the beginning of each run allows for 
introducing only the frequencies of interest during each 
part of the mission. 



Fig. 2 
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When, during deployment, the need for considering 
another frequency appears, the simulation requires a restart 


with a new 

description of the 

tether 

which 

takes 

into 

account new 

mass points. 

The 

I.C. 

for 

such a 

new 

conf iguration 

can be derived 

f rom 

the 

results of 

the 


previous simulation, under the assumption that the new modes 
are not excited; thus, the new masses must be added until 
the associated mode effects are negligible. Reverse 
considerations can be applied for the retrieval phase. 

An alternative/complementary model is being developed, 
to be used for cross validation purposes: the mass points 
are maintained uniformely spaced along the tether length to 
represent the actual mass distribution, during both 
deployment and retrieval. The system frequencies are now 
closer to the actual ones, but the mass points "move" along 
the tether. 

2.3 System Model 

The Orbiter and the Satellite are considered rigid 
bodies connected by a tether represented by series of 
springs and point masses. The system is considered orbiting 
around the Earth. 

The references frames, showed in figure 3, are: 

- Inertial reference frame (IRF), located in the centre 
of the Earth. 
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- Accelerating reference frame <.ARF), in the centre of 
mass of the overall system; it is always parallel to the 
IRF. The motion of ARF origin is described with respect to 
IRF, integrating the orbit equations of the system. The 
motion of the parts of the system is described with respect 
to ARF. This is to avoid possible numerical problems when 
the motion of a system is referred to a very far reference 
frame (small numbers added to large numbers). 

- Orbiter body reference frame, located in the Orbiter 

c . m. 

- Satellite body reference frame, in the Satellite 
c.m.; z-axis i3 directed toward the tether attachment point. 

- Satellite attitude reference frame, in the Satellite 
c.m., with the z-axis always directed toward the Earth 
centre. The Satellite attitude is referred to this frame 
( local vertical). 

2.4 Dynamic Equations 

The motion of each body is computed with respect to 
ARF. The Orbiter has 3 translational rigid d.o.f's; the 3 
rotational d.o.f.'s can be taken into account or neglected, 
according to the figure of an input flag. The Satellite has 
all the 6 rigid d.o.f.'s. The mass points describing the 
tether have only 3 translational d.o.f's. 
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Three differential equations describe the orbital mo- 
tion of the origin of MF about the Earth. 

The translational differential equations with respect 
to ARF are obtained for each body considering all the forces 
acting on it (actual and apparent ones): 

mx = -mA + mg + T + F 

where: X : body c.m. position in ARF 

~k : acceleration of ARF origin in IRF 
g : local gravity acceleration 
T : force(s) coining from the tether 

F : other external forces (thrusters firing, imposed 
perturbations) . 

No Coriolis force is present because the ARF does not 
rotate. The rotational differential equations for the 
Orbiter or the Satellite are represented by the Euler equa- 
tions. The rotational motion is described by means of 3 
Euler angles, deriving from the integration of the Euler 
equations (motion about the centre of mass). 

2.5 Other Features 

No environmental model has been implemented: aerodynamic 
drag and electromagnetic forces on the tether were consi- 
dered negligible, even if the possibility to input forces on 
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the mass points is foreseen. 

In general arbitrary forces can be introduced to 
analyse the system behaviour under specified perturbations. 

All the Satellite thrusters are modelled and their 
misalignements , as source of attitude perturbations, consi- 
dered. 

The control laws of the tether length and in-plane & 
out-of -plane libration developed for the TSS program have 
been implemented. 

The Orbiter attitude may have an indirect but important 
impact on the Satellite behaviour during some mission 
phases, so it has been included in the model together with a 
very simplified model of the digital autopilot and reaction 
control system. 

Due to lack of data on tether characteristics, no 
damping has been considered, for the system oscillation, this 
resulting in a conservative hypothesis. The tether torsional 
stiffness is negligible and so it is not taken into account 
in this program. 

A new feature which will be introduced in the program 
is represented by the addition of flexible modes on the 
Satellite , simulating two deployable booms. 

The need to model in an easy way the flexibilities of 
the deployer boom on the Orbiter, expecially at very short 
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tether length. Is being investigated. 

3 - SOFTWARE DESCRIPTION 

The general architecture of the program is presented in 
figure 4. The simulation program is coded in FORTRAN H and 
runs on a VAX 785 machine. 


Fig. 4 
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The YDT routine computes, time by time, the state 
vector derivative to be integrated. It determines all the 
force and torque vectors acting on each part of the system, 
via the management of the call of several subroutines, each 
describing a particular item (see table 1), allowing to build 
the differential equations. 

Subroutine Function 

LOCVERT Sat. attitude wrt local vertical computation 

AMCS AMCS model 

ACS Simplified AMCS model 

INLINE In-line thrusters command logic 

PTCHCMD Commanded libration pitch angle 

LIBRCTL Libration control logic 

ORBATT Orbiter attitude control system model 

TENSCMD Commanded tether tension 

REELEQ Tether reeling equation 

VARMAS Tether variable characteristics computation 

THRUST Thrusters forces and torques computation 
ORBEQS Orbital equations 

TENSION Tether tension force acting on each body 
GRAVITY Gravity forces computation 

TRANEQS Translational equations for each body 
ROTEQS Rotational equations (Satellite , Orbiter ) 

TABLE 1 

c-^ 
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The program has been conceived in a modular way so that 
the change of the generic module can be friendly performed, 
as well as new routines introduced and other bypassed. 

The integrator is a 4th order Runge-Kutta algorithm. 

The output of the program are: 

- diagrams of selected variables; 

- print-outs; 

restart file, which allows for the restart of the 
simulation from the end of a previous one; 

- output file , stored on tape, containing the time 
history of all the variables needed as input for the 
verification and testing programs. 

The dynamic part of the program is being checked via 
comparison runs to be performed with a qualified simulation 
program. 
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INTERMEDIATE STOPS DURING DEPLOYMENT AND/OR RETRIEVAL 


TETHERED SATELLITE SYSTEM FREQUENCY SPECTRUM 
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PARTICIPATION IN PAYLOAD INTEGRATION PLAN (PIP) ANNEX TWO SUPPORT 
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TETHER FRAME (ORIGIN AT ONE ATTACH POINT, Z-AXIS ALONG TETHER LENGTH) 
NADIR FRAME (ORIGIN AT CENTER-OF-MASS, Z-AXIS ALONG RADIUS VECTOR) 
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19 

0.0928 

9879E+02 

• 

20 

0.0977 

8472E+02 


21 

0, 1025 

-. 8446E+02 


22 

0. 1074 

-. 8436E+02 


23 

0.1123 

-. 7976E+02 


24 

0. 1172 

-.7733E+02 

X 

23 

0.1221 

861 2E+02 

26 

0.1270 

-. 8866E+02 


27 

0. 1318 

-, 8566E+02 

S> 

28 

0. 1367 

-. 1 00CE+03 

29 

0. 1416 

- . 874 1 E+02 


30 

0. 1465 

- , 866 1 E+02 

} 

31 

0. 1514 

-. 1 000E+03 

32 

0. 1563 

-. 81 08E+02 


33 

0. 1611 

-.3741E+02 

D 

34 

0. 1660 

- . 382 1 E + 02 

33 

0.1709 

-. 370 1E+02 - 


36 

'571755* 

5449E *02 


37 

0. 1807 

-.71 46E+02 

• 

38 

0. 1855 

-. 8704E4-02 

1 

39 

0. 1904 

-. B867E+02 

<S 

40 

0. 1953 

-. 9200E+02 

41 

0.2002 

-. 1000E+03 


42 

0.2051 

-.9790E+02 


43 

0.2100 

-. 8 1 43E+02 

J 

44 

0.2148 

-. 8039E+02 


43 

0.2197 

-. 8832E+02 


46 

0.2246 

-. 9532E+02 


47 

0.2295 

-. 7533E+02 


48 

0.-2344 

-.7613E+02 

. 

49 

0.2393 

-. 8338E+02 

N-/ 

30 

0.2441 

-. 9637E+02 


51 

0.2490 

-. 1000E+03 


52 

0.2339 

-. 9682E+02 


53 

0.2588 

-. 9738E+02 


54 

0.2637 

8560E+02 

J 

55 

0.2686 

-. 81 48E+02 

36 

0. 2734 

-.0144E*O2 


37 

0.2783 

B889E+02 
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j 
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55 

0.2686 

56 

0.2734 


57 

0.27B3 


58 

0.2832 

• 

59 

0.2881 


60 

0.2930 

* 

61 

0.2979 

62 

0.3027 


63 

0.3076 


64 

0.3125 

* 

65 

0.3174 


66 

0.3223 


67 

0.3271 

• 

68 

0.3320 


69 

0.3369 

) 

70 

0.3418 

71 

0.3467 


72 

0.3516 


73 

0.3564 

• 

74 

0.3613 


75 

0.3662 


76 

0.3711 


77 

0.3760 


78 

0.3809 


79 

0.3857 

• 

80 

0.3906 


81 

0.3955 


82 

0.4004 


83 

0.4053 


84 

0.4102 


85 

0.4150 


86 

0.4199 


87 

0.4248 


88 

0.4297 

S 

89 

0.4346 


90 

0.4395 


91 

0.4443 

• 

92 

0.4492 


93 

0.4541 


94 

0.4590 

95 

0.4639 


96 

0.4680 


97 

0.4736 

• 

98 

0.4785 


99 

0.4834 

*, 

100 

0.4883 

• 

101 

0.4932 


102 

0.4980 


103 

0.5029 

• 

104 

0.5078 


105 

0.5127 


106 

0.5176 

107 

0.5225 


108 

0.5273 


109 

0/5322 


110 

0.5371 


111 

0.5420 


112 

0.5469 

• 

113 

0.5518 


114 

0.5566 


115 

0.5615 

A 

116 

0.5664 


417 

0.5713 


118 

0.5762 


r.O0CH't*U* 

-.81 48E+02 KAdtAU C Or<r. 

- . 8 1 4 4E + 02 
“. 8889E+02 
-. 8990E+02 
8865E+02 
-.9329E+02 
- . 9269E *02 
- . B059E+02 
-.10006+03 
8425E +02 
-.89906+02 
-.93566+02 
-. 85146+02 - J 
-.89276+02 
-. 1000E + 03 
8535E+02 
-.75446*02 
- . 8493E+02 
- . 8754E+02 
8942E+02 
1 000E+03 
-.881 OE + 02 
-. S925E+02 
-. 9559E+02 
- . 8688E+02 
-. 8564E+02 
1 0006+03 
8257E+02 
-. 8303E+02 
1 OOOE + 03 
- . B677E+02 
- . 8600E+02 
10006+03 
8373E+02 
- . B266E+02 
1000E+03 
-. 1 000E+03 
1000E+03 
- . 870 1 E *02 
- . 8441 E*02 
-. 1 OOOE*C3 
-.75046+02 
-. 736 L E*02 
- . 9064E+02 
- . 9 1 56E+02 
8756E+02 
9458E+02 
-.931 0E+02 
10006 + 03 
1 000E+O3 
8436E+02 
8393E+02 
-.92586+02 
-, 7683E+02 
-. 8798E+02 
-.93426+02 
-. 79B4E+02 
-. 8506E+02 
-.90316+02 
-.92206+02 
-.10006+03 
9986E+02 
-.83716+02 
-.861 4E+02 
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***** LOG. POWER SPECTRUM 


<2 = i.El fi. 


LIN. SPETTRALE FREQUENZA 


VALORE (DB) 


( £ v£iO fj 


• 

.0 

0.0000 

0.26406*03 


1 

0.0049 

0. 2520E+03 

tt 

2 

0.0098 

0.11 08E*03 


3 

0.0146 

0. 1276E+03 

* 

4 

0.0195 

0. 15306*03 

3 

0.0244 

0. 1988E*03 


6 

0.0293 

0. 20336*03 


7 

0.0342 

0. 1871E*03 

• 

B 

0.0391 

0. 1406E*03 


9 

0.0439 

0. 1190E*03 


10 

0.0488 

0.10406*03 


11 

0.0537 

0. 92316*02 


12 

0.0586 

0. 82476*02 


13 

0.0635 

0 . 7373E*02 

Jt 

14 

0.0684 

0.6548E+02 


15 

0.0732 

0. 5724E*02 


16 

0.0781 

0. 4852E*02 

• 

17 

0.0830 

0.41 40E*02 


IB 

0.0879 

0. 42936*02 


19 

0.0928 

0. 49326*02 

• 

20 

0.0977 

0.35476*02 


21 

0. 1025 

0. 6098E+02 


22 

0.1074 

0.6615E+02 

• 

23 

0.1123 

0. 7 l 24E*02 


24 

0.1172 

0. 76476*02 


23 

0.1221 

0.8201E*02 

• 

26 

0.1270 

0. 88016*02 


27 

0. 1318 

0 .94676*02 


28 

0.1367 

0.10226*03 


29 

0. 1416 

0. 1 1 10E*03 


30 

0. 1465 

0.121 6E*03 


31 

0. 1514 

0. 13516*03 

• 

32 

0. 1563 

0. 1536E+03 


33 

0. 1611 

0. 1852E + 03 


34 

0. 1660 

0.21646*03 

33 

0.1709 

0.21 826*03 


36 

0. 1758 

0. 19346*03 


37 

0.1807 

0. 1572E+03 

• 

38 

0. 1855 

0. 1 375E + 03 


39 

0. 19*4 

0. 1235E + 03 


40 

0./V33 

0.1126E*03 

• 

41 

Q. 2002 

0. 1 033E*03 


42 

U. 2051 

0. 9381E+02 


43 

0.2100 

0 . 8904E*02 

• 

44 

0.2148 

0.8299E*02 


43 

0.2197 

0. 7746E*02 

A 

46 

0.2246 

0. 7230E*02 

9 

47 

0.2295 

0 , 6737E+02 


48 

0.2344 

0. 62536*02 


49 

0.2393 

0. 37616*02 

# 

30 

0.2441 

0.52326*02 


31 

0.2490 

0. 46406*02 


32 

0.2539 

0. 39498*02 

• 

53 

0.2588 

0.34018*02 


34 

0.2637 

0. 37376*02 

•> 

33 

0.2686 

0.43686*02 

56 

0.2734 

0. 33336*02 


57 

0.2783 

0.6027E*02 


58 

0.2832 

0.66966*02 
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av 

60 

61 

62 

63 

64 

65 

66 
67 

<j. /aw i 
0. 2930 
0.2979 
0.3027 
0.3076 
0.3125 
0.3174 
0.3223 
_0. 3271 

0 . 

0.81 05E *02 
0.891 JE*02 
0 . 9844 E + 02 
0. 1097E+03 
0. 1 240E+03 
0. 1 442E+03 
0. 18I0E+03 
0. 2023E+03 

68 

0.3320 

0. 1987E+03 

6? 

0.3369 

0. 1658E+03 

70 

0.3418 

0. 1369E+03 

71 

0.3467 

0. U90E + 03 

72 

0.3516 

0. 1059E+03 

73 

0.3564 

0. 9349E + 02 

74 

0.3613 

0. 8673E+02 

75 

0.3662 

0. 79I6E+02 

76 

0.3711 

0. 7235E+02 

77 

0.3760 

0.6606E+02 

78 

0.3809 

0. 6001E+02 

79 

0.3857 

0.3390E+02 

80 

0.3906 

0. 4732E+02 

81 

0.3955 

0. 3940E+02 

82 

0.4004 

0.2796E+02 

85 

0.4053 

0. 6709E+01 

84 

0.4102 

0. 2424E+02 

85 

0.4150 

0. 36B8E+02 

86 

0.4199 

0. 4503E+02 

87 

0.4248 

0. 5149E+02 

88 

0.4297 

0. 3729E+02 

89 

0.4346 

0. 6287E+02 

90 

0.4395 

0. 6B33E+02 

91 

0.4443 

0.7449E+02 

92 

0.4492 

0. 8095E+02 

93 

0.4541 

0. 8813E+02 

94 

0.4590 

0. 963SE+02 

95 

0.4639 

0. 1061E+03 

96 

0.4688 

0. 1180E+03 

97 

0.4736 

0. 1336E+03 

98 

0.4785 

0. 1562E+03 

99 

0.4834 

0. 1886E+03 

100 

0.4883 

0. 1991E+03 

101 

0.4932 

0. 1880E+03 

102 

0.4980 

0. 1535E+03 

103 

0.5029 

0. 1332E+03 

104 

0.5078 

0. 1 177E+03 

105 

0.5127 

0. 1058E+03 

106 

0.5176 

0. 9606E+02 

107 

0.5225 

0. 8777E+02 

108 

0.5273 

0. 8047E+02 

109 

0.5322 

0.7381E+02 

110 

0.5371 

0. 6754E+02 

111 

0.5420 

0. 61 38E+02 

112 

0.5469 

0. 5494E+02 

113 

0.5518 

0. 4773E+02 

114 

0.5566 

0. 3891E+02 

115 

0.5615 

0. 2933E+02 

116 

0.5664 

0. 3302E+02 

117 

0.5713 

0. 4243E+02 

118 

0.5762 

0. 3004E+02 

119 

0.5011 

0. 363IE+02 

120 

0.5859 

0. 6207E+02 

121 

0.5908 

0.6763E+02 

122 

0.5957 

0. 7 J23E+02 

123 

0.6006 

0.7916E+02 


119 

0.5811 

0.563SE+02 

120 

0.5859 

0 .6207E+02 

121 

0.5908 

0.6763E+02 

122 

0.5957 

0. 7323E+02 

123 

0.6006 

0 . 791 6E+02 

124 

0.6055 

0. 8556E+02 

123 

0.6104 

0.9266E+02 

126 

0.6152 

0. 100BE + 03 

127 

0.6201 

0. 11Q2E*03 

128 

0.6250 

0. 1217E + 03 

129 

0.6299 

0. 1364E + 03 

130 

0.6348 

0.1561 E+03 

131 

0.6396 

0. 1 795E + 03 

132 

-0.6445 

0. 1892E+03 

133 

0.6494 

0. 1810E+03 

134 

0.6543 

0. 1 58 IE *03 

135 

0.6592 

0. 1 378E+03 

136 

0.6641 

0. 1 229E+03 

137 

0.6689 

0.11 12E+03 

138 

0.6738 

0. 1016E+03 

139 

0.6787 

0 . 9356E+02 

140 

0.6836 

0. 8655E+02 

141 

0.6885 

0 . 8036E+02 

142 

0.6934 

0.7480E+02 

143 

0.6982 

0. 6975E+02 

144 

0.7031 

0. 6506E+02 

145 

0.7080 

0. 6069E+02 

146 

0.7129 

0.5651E+02 

147 

0.7178 

0. 5248E+02 

148 

0.7227 

0. 4844E+02 

149 

0.7275 

0. 4428E+02 

150 

0.7324 

0.3976E*02 

151 

0.7373 

0. 3442E+02 

132 

0.7422 

0. 2755E+02 

153 

0.7471 

0. 1916E + 02 

134 

0. 7520 

0. 2098E+02 

155 

0.7568 

0.3318E+02 

156 

0.7617 

0. 4394E*02 

157 

0.7666 

0. 5376E+Q2 

158 

0.7715 

0 . 6377E+02 

139 

0.7764 

0 . 7490E+02 

160 

0.7813 

0. 8805E+02 

161 

0.7861 

0. 1041E + 03 

162 

0.7910 

0. 1 2 1 5E+03 

163 

-CL. 7959 

0. 1311E+03 

164 

0.8008 

0. 1275E + 03 

165 

0.8057 

0. 1124E+03 

166 

0.0105 

0. 951 5E+02 

167 

0.8154 

0. 8079E+02 

168 

0.8203 

0. 6903E+02 

169 

0.8252 

0.5890E+02 

170 

0.8301 

0.4968E+02 

171 

0.8350 

0. 4079E+02 

172 

0.8398 

0. 321 4E+02 

173 

0.8447 

0. 2645E+02 

174 

0.8496 

0. 2815E+02 

175 

0.8545 

0. 3312E+02 

176 

0.8594 

0. 3821E+02 

177 

0.8643 

0.4277E+02 

178 

0.8691 

0. 4700E+02 

179 

0.8740 

0.51 19E+02 

180 

0.8789 

0.5535E+02 

181 

0.8838 

0.5964E+02 

182 

0.8887 

0. 64 18E+02 

183 

0.8936 

0.6898E+02 
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182 

0.8B97 

0 . 64 l 8E+02 


183 

0,8936 

0. 6899E+02 


184 

0,8984 

0. 7424E+02 


185 

0,9033 

0. 79996+02 


186 

0 , 9082 

0.8639E+02 


187 

0.9131 

0. 9363E+02 

188 

0.9180 

0. 1019E+03 


189 

0.9229 

0. 11166*03 


190 

0.9277 

0.12296*03 

J 

191 

0.9326 

0. 1360E+03 


192 

0.9375 

0.14896*03 


193 

0.9424 

0. 1363E+03 

194 

“6.9473 

~ 6. 1548E*03 _ 


195 

0.9521 

0.14476*03 


196 

0.9570 

0. 131 3E*03 

197 

0.9619 

0. 1188E+03 


198 

0.9668 

0. 1080E*03 


199 

0.9717 

0.98786+02 

200 

0.9766 

0. 9075E+02 


201 

0.9814 

0.83626*02 


202 

0.9863 

0.771 7E+02 

203 

0.9912 

0. 7 1 1SE+02 


204 

0.9961 

0. 6543E+02 


203 

1.0010 

0. 600 IE+02 

206 

1.0059 

0. 55 1 2E+02 


207 

1.0107 

0. 51 B8E+02 


208 

1.0156 

0.5196E+02 

209 

1.0205 

0. 5S20E+02 


210 

1.0254 

0. 5993E+02 


211 

1.0303 

0.6513E+02 

212 

1.0352 

0. 705 1E+02 


213 

1.0400 

0. 761 7E+02 


214 

1.0449 

0. 821 7E+02 

215 

1.0498 

0. 886SE+02 


216 

1.0547 

0. 9S89E+02 


217 

1.0596 

0. 1 040E+03 

218 

1.0645 

0.11 30E+03 


219 

1.0693 

0. 1232E+03 

% 

220 

1.0742 

0. 1 337E+03 

221 

1.0791 

0. 1426E+03 


222 

1.0840 _ 

_ 0.1466E+03 


223 

1.0889 

0. 1440E + 03 

224 

1.0938 

0. 1 358E+03 


225 

1.0986 

0. 1254E+03 


226 

1.1035 

0. 1 151E+03 

227 

1.1084 

0. 1058E+03 


228 

1.1133 

0. 9765E+02 


229 

1.1182 

0. 9041E+02 

230 

1.1230 

0.B401E+02 


231 

1.1279 

0. 7831E+02 


232 

1.1328 

0.731 6E+02 

233 

1.1377 

0 . 6856E+02 


234 

1.1426 

0. 64 46E + 02 


235 

1.1475 

0. 6083E+02 

236 

1.1523 

0. 5777E+02 


237 

1.1572 

0. 5543E+02 


238 

1.1621 

0. 5402E+02 

239 

1.1670 

0.5396E+02 


240 

1.1719 

0. 5563E+02 

% 

241 

1.1768 

0. 5897E+02 

242 

1.1816 

0, 6386E+02 


243 

1.1865 

0. 6997E+02 


244 

1.1914 

0. 77 10E+02 


245 

1.1963 

0. 8508E+02 
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247 

1.2061 

0. 1 023E *03 

248 

1.2109 

0. 1097E+03 

249 

1.2158 

0. 1I3SE + 03 

250 

1.2207 

0. 1 1 32E + 03 

251 

1.2256 

0. 1081E+03 

252 

1.2303 

0. 1003E+03 

253 

1.2354 

0.9160E+02 

254 

1.2402 

0.831 4 E ♦ 0 2 

255 

1.2451 

0.7541E+02 

256 

1 . 2500 

0. 6860E *02 

257 

1.2549 

0 . 6293E +02 

258 

1.2598 

0. 5860E + 02 

259 

1.2646 

0.5594E+02 

260 

1.2695 

0.5505E+02 

261 

1.2744 

0. 5582E+02 

262 

1.2793 

0.5782E+02 

263 

1.2842 

0.6071E+02 

264 

1.2B91 

0. 6428E+02 

265 

1.2939 

0.6839E+02 

266 

1.2988 

0. 7301E+02 

267 

1.3037 

0. 7812E+02 

268 

1.3086 

0. 8373E + 02 

269 

1.3135 

0 . 8987E+-02 

270 

1.3164 

0.9651E+02 

271 

1.3232 

0. 1035E + 03 

272 

1.3281 

0. 1106E + 03 

273 

1.3330 

0. 1170E+03 

274 

1.3379 

0. 1216E+03 

275 

1.3428 - 

_ _0.1232Eo.Ql 

276 

1.3477 

0. I214E+03 

277 

1.3525 

0. 1 166E+03 

278 

1.3374 

0. 1102E+03 

279 

1.3623 

0. 1032E+03 

280 

1.3672 

0.9633E+02 

281 

1.3721 

0 . 8990E + 02 

282 

1.3770 

0.8400E+02 

283 

1.3818 

0.7874E+02 

284 

1.3867 

0.741 7E + 02 

285 

1.3916 

0.7045E+02 

286 

1.3965 

0.6789E+02 

287 

1.4014 

0.6691E+02 

288 

1.4063 

0. 6740E+02 

289 

1.4111 

0.6954E+02 

290 

1.4160 

0 . 7290E + 02 

291 

1.4209 

0.771 2E+02 

292 

1.4258 

0 . 8200E+02 

293 

1.4307 

0. 8739E+02 

294 

1.4355 

0. 9315E+02 

295 

1.4404 

0. 9902E+02 

296 

1.4453 

0 . 1 046E+03 

297 

1.4502 

0. 1092E+03 

298 

1.4351 

0.111 9E+03 

299 

1.4600 

0.11 22E+03 

300 

1.4648 

0. 1100E + 03 

301 

1.4697 

0. 1057E+03 

302 

1.4746 

0. 1001E + 03 

303 

1.4793 

0. 9399E+02 

304 

1.4844 

0. 0772E+O2 

305 

1.4893 
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0 

0,0000 

0. 57 40E-01 


1 

0.0049 

0. 2439E-01 


2 
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a 

3 
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4 
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5 

0.0244 
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• 

6 
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7 
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• 

9 

0.0439 

0. 1 1 32E+00 


10 
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n 
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\ 

12 
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13 
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15 
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17 
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• 

18 
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19 
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20 
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• 

21 
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22 
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0. 1 531E-01 


23 

0.1123 

0. 1 574E-01 
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24 

0.1172 

0. 1 037E-01 


25 
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0. 1 278E-01 


26 
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27 
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28 

0.1367 
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29 
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30 
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31 
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0.4397E-02 


32 
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33 
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34 
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35 
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44 
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• 

45 
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46 
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0, 201 6E-01 


47 
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• 

48 
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0, 1643E-01 


49 
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0. 4621E-02 


50 

0.2441 

0. 1 060E-01 

• 

51 
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0. 7 1 32E-02 


52 
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53 
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0.8347E-02 

• 

54 
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55 
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0. 151 0E- 01 


56 
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57 
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58 
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0.821 6E-02 


63 
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0.22118-01 

64 
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65 

0.3174 
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66 
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67 
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£8 
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69 
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70 
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71 
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73 
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74 
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75 
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76 

0.3711 
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77 
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78 
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79 
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60 
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82 
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83 

0.4053 

0. 1246E-01 


84 
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85 
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87 
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90 
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91 
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92 
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93 
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94 
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95 
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96 
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97 
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98 
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0.5518 
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0.5566 
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115 

0.5615 
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116 

0.5664 
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117 

0.5713 

0.11 45E-01 

118 

0.5762 

0. 1556E-01 


119 

0.5811 

0. 6703E-O2 


120 

0.5059 

0. 9616E-02 

121 

0.5908 

0.9919E-02 


122 

0.5957 

0. 3327E-02 


123 

0.6006 

0.9681E-02 


124 

0.6055 

0. 6503E-03 


125 

0.6104 

0. 1298E-01 
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125 
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127 
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128 
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• 
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0.1391E-01 
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0.4410E-02 


134 

0.6543 

0. 1043E-01 
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9 
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140 

0.6836 
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142 
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143 

0.6982 
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144 

0.7031 

0. 1944E-01 


145 

0.7080 
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146 
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147 

0.7178 
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148 
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149 
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0. 4 1 6 IE-02 
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0.7262E-02 
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159 
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160 
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0.7861 
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164 

0.8008 
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0. 1 160E-01 


166 
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0. 1460E-01 


167 
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169 
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0.8390 

0.1484E-01 


173 

0.0447 

0.2574E-01 
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0.8496 

0. 1780E-01 
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ABSTRACT 


Dynamics of Tethers in a linearized analysis can be 
considered as the superposition of propagating waves. 
This approach permits to have a new way for the analysis 
of a Tether behaviour during deployment and retrival, 
were a Tether can be considered composed by a part at 
rest and a part sujected to propagation fenomena, being 
the separating section depending on time. 

The dependence on time of the separating section 
requires the analysis of the reflection of the waves 
travelling toward the part at rest. Such a reflection 
generates a reflected wave, whose characteristics are 
determined. 

The propagation fenomena of major intest in a Tether are 
transverse waves and longitudinal waves, all 
mathematically modelled by the "vibrating chord" 
equations, if the tension is considered constant along 
the Tether itself. An interesting problem also 
considered is concerned with the dependence of the 
Tether tension from the longitudinal position, due to 
microgravity, and the influence of this dependence on 
propagation waves. 
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INTRODUCTION 


Dynamics of Tethers(')* as well as of any structure, in a linearized 
analysis can be considered as the superposition of propagating waves. 
This requires the study of dynamic propagation along the Tether with 
the appropriate boundary conditions. During deployment and retrieval, 
with reference to a lagrangian reference system ("), Tether can be 
considered as composed of two parts - one at the rest and one subjected 
to the propagation phenomena. These two parts are separated by a 
section that changes with time, i.e. the Tether section that bounds the 
part constrained to rest is changing with time (the other part being 
free to move and vibrate and having the opposite end section subjected 
to the boundary conditions imposed by the satellite). 

The propagation phenomena of major interest from a practical point of 
view are the following. 

1) Transverse waves, mainly a "vibrating chord" behaviour, where 

inertial forces and the tension in the Tether - in combination with 
its local curvature - are the most important elements of the 
dynamic equilibrium. 

2) Longitudinal waves, mathematically modelled by the "vibrating chord" 
equation, where inertial forces and longitudinal internal forces, 
due to elastic deformations, are the most important elements of the 
dynamic equilibrium. 

The dependence on time of the section which bounds the part at rest 
requires the analysis of the reflection of the wave travelling toward 
the part at rest. Such a reflection generates a wave travelling 
outward, whose characteristics are to be determined. 

The A. had previously considered from a theoretical point of view such 
problem in particular in order to analyze the behaviour of deployable 
booms subjected to longitudinal and flexural dynamic phenomena (see 
(1), (2), (3) and (4)). Also in the case of the problems concerning a 

structure like the Tether, the theoretical analysis gives rigorous 
solutions and permit an insight into experimentally observed effects, 


(') It can be suggested Ref. (8) for a general presentation of Tether 
concept and its (dynamic) problems. 

(") A reference system which introduces a bi-univocal correspondence 
between a longitudinal coordinate and each Tether section. 
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that by some authors were erroneously thought to be '’continuous" 
changes of frequencies and amplitudes of the proper modes. 

This paper belongs to a series of works having the scope of opening a 
new way in the approach by means of mathematical models of several 
mechanical problems of deployable systems of telescopic and Tether 
type. 

Recently several attempts have been made to solve the problem of the 
telescopic structures behaviour. Such attempts are mainly based in 
changes of the coordinates in order to take into account the changes 
during the time of the space where the problem is defined. 

As a matter of fact these attempts don’t seem obtain good results. They 
don't take into account energetic balances. 

On the contrary this work introduces and developes to some extent the 
basic idea of considering each dynamic motion in a structure as the 
results of wave propagation, taking also into account energetic 
exchanges at the ends. 

In the case of vibrating chord the problem of the time dependence of 
the definition space can't be resolved by means of the Cauchy, Goursat 

and Darboux results, (5). These results deal with the problems of time 

depending location of the sections where are imposed the boundary 
conditions. In the part that is external to such sections dynamic 

phenomena take place that are coherent with phenomena acting in the 
internal part and contribute to supply or spillover energy in it. 

This work deals with the request of having external parts at the rest 
(not only the boundary sections). Therefore boundary sections have the 
behaviour of surface where internal dynamic phenomena "reflect". 
Obviously reference is made to a constant section unaxial structure, as 
Tether can be considered. Longitudinal tension loads due to 

microgravity permits to consider additive small tension or compression 
loads without critical phenomena. 
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FUNDAMENTAL EQUATIONS 


Let us have first a brief recall of concepts, with reference only to 
the case of longitudinal waves making use of Ref. (1). The reader can 
easily do the extension to transverse waves. The problem can be 
analyzed by means of an equilibrium and a continuity equations. 

If 6-^X/t) is the stress at point X and time t and M. is 
the velocity of the motion, the equilibrium linearized equation is 


1) 


9 ) 6 " 

- - P 

9t -5 9t 


The continuity linearized equation on the other hand is 

2) A Qg 

’dx ~ E Or 


If we put c ) 


l’) 


2 ' ) 


operating we obtain 

A 0 2 t 


3) 


4) 


Car 

*• c 

; z 
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± 


9^ 

2 * « 

Z 2 ' 

W 2 - 

solutions of eq.i 

s 1* 

) and 

Gr = 

h (*- 

<L 

J ^ 

M* - 

^A* 

X 

C 

J * 


Eqs. 3) and 4) indicate that the motion of the bar is composed of two 
waves: one travelling in the increasing x direction and another 
travelling in the decreasing x direction. 

On the base of these considerations it is possible to obtain a relation 
between p and u of each travelling wave. 
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t from 3) and 4) one obtains 


If 4. = i- r and 


a 


5) 


0 ^ 

£>r 


^ 

c. 

, jf. I 

C ■ Otv J 


S -JM. 

*i 

9t# 
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V ' '"2- 

Taking into account eq. 1) from 5) and 6) we have the relation 
7) 


m-u ci^l;uuuu cq. i / u um u / euiu u; we nave me 

*/ / ^A/ 1 _ / o/ \ 

T( 77 - IT. * > 77 77: / ' 


**??, y h 

If A , eq. 7), gives the relation 

'* ' * A _ pc **- _ £. |i» 

*», ~ J At, ^ 3?< 

On the other hand if j\, « t> , we have the relation 

dfii _ 4 <Lji z _ £_ 

itl J Zf, c 


With the initial conditions /t'j J ^ ^ ) - O or (separately) 

[io)=OM 2 (t^)=Othe following relations between u and p hold for each 
travelling wave 


^(*'71 * | 


Such relations enable us to determine one of the two values ^ or when 
the other is known (eq. 8) and also to determine one of the values or 
when the other is know (eq. 9). 

They are connected respectively with the separate behaviour of the two 
travelling waves. 
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REFLECTION ON THE TIME DEPENDENT SECTION 


Also in this paragraph longitudinal waves are considered as a sample 
problem, still making use of Ref (1). The extension to transverse wave is 
easy. In practical applications the boundary conditions that are usually 
considered for uniaxial extensional bars are free edge (<T-^ ) an d fixed 

edge (-4*0). At the edge where deployment is done, it is (OOc c )=0, where 
) i s the time dependent section that can be considered as fixed 
and =? \T is the velocity of displacement of the constraint. Here 

is an abscissa on the indeformed bar. 

As a first analysis of the behaviour at a time dependent fixed section we 
can consider the problem of an extensional bar having a free edge at 
and constrained with M/tx^) =• O = d 0 -t-CV % where l , <T are 

constant. The bar is subjected in X-O to an external extensional 
specific force indipendent of the timet. 

Such force produces a wave travelling in the direction of the increasing 
X when such wave reaches the moving constraint section a 

reflected wave of specific force is generated that runs in the 

direction of the decreasing x. 

We will now determine the characteristic of the reflected wave, before it 
reaches the section X=0 

In order to determine the reflected wave by means of an energy balance is 
necessary to dispose of an evaluation of the energy exchanged at the 
constrained end, A discussion on this subject is performed in ref. (1). 

The conclusion is that the constraint has and energy exchange different 
from zero e that the reflected wave can be determined by means of a 
behavioural analysis like the following. 

During the time we have the already introduced displacement of the 
section where the constraint is imposed. Such displacement corresponds to 
the internal deformation of the rod. 

The interval d.Tc^'fdX during df , withstands a length change dt = \rft~ 

with a strain — 4,1 , ‘ ' 

On the other hand the final stress in rctt 
= 6 t r < ^) — , that means a strain £ s C-&, 


must be 5" - ^ + fl^ = 
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Equating the two expressions of the strain gives 


10 ) 


- 




C - 

C * xT 


This result 10) coincides with that of the application of the Goursat 
Darboux and Cauchy problem solution, ( 5) , with the condition atx- 
In spite of the observed coincidence with the well known results of the 
vibrating chord analysis, the proposed model presents the advantage of 
the applicability to more complex problems as dispersive systems, (see 
for instance (6) and (3)). 
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EXTENSION TO TIME DEPENDENT AMPLITUDE (OF THE TRAVELLING WAVE) AND SPEED 
(OF THE BOUNDARY CONDITION) 


Let us consider in a non dispersive system the wave travelling 

inward the time variable restraint and the travelling outward 

expanded as follows, (see also (2)), 

+ 0(,(TC-Ct) 4 . ... , 

JA Z * -t /3 (x+ci J - - ' • 

If at a time t the restraint condition is 

2 i«o| at the time ~t + At where At will tends to zero, 

without loss in generality we can put O and introduce a dummy 

variable t such that £ < t J-At. 

During AtT at the section the displacement 

12) Ai *Jj\ (D cj 4 . uc 7 (o, t)]it 

take place. During the same in the region X>0 the internal 

deformation generates a change in length 

13) A i (JZ; At ) _ it. 

The equation between such two i substituting 14) and taking into 

account only the terms of the lower order for respect to the principal 
At , operating gives^ 

C^z) A> = J 

or 

Relation 14) means that, beside the higher order terms, 
like in the case of and-tt^of constant value. 
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EXTENSION TO PROPAGATION SPEED DEPENDING ON THE LONGITUDINAL POSITION 


The present paragraph concerns the extension of the analysis to the case 
where eq. 1’) and 2') become as follows 


9% 

A 

9*? 

A") 

TP ~ 

c 7 fr) 

9 1*- ; 
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%‘j 

TP 

£ Z (X) 

Of 2 - * 


This case includes the dynamics of transverse waves in Tethers where 
longitudinal load is depending on the position due to microgravity. 

Because ct'xj is not constant but a function of oc / C -<&)> expressions 
3) and 4) are not still valid, and it is necessary to find an appropriate 
way of solution. 

The A. in a previous paper, (7) here largely recalled, proposed that the 
general solution of eq.s 1") were composed by means of two waves in 
opposite directions travelling and having speed depending on x. Such 
waves reduce to eq.s 3) and 4) when reduce to a constant C. If<^ 

indicates or -to as necessary, the following expression was adopted (') 


± 4 cix) 


AS) 


could be the 


where is an arbitrary constraint and 
propagation speed at ^ . 

Obviously, not all the functions and <j)^ are usefull to satisfy 1"). The 
problem is now reduced to the determination of the functions and <j^if 
any, that can satisfy 1" ). 

The functions A and b can be examined separately. 

^ = t-J hf> • we have 

• (4i 9i dd 4 / _ i 

^ - ku*)J 


If we put 

9* 2 
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ITT di Q t di 


1 14 i (j__ 

k c l fr) 
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d i 


( * ) When ^^Jsconst. and k-'H , the proposed expression reduces to 3) 
and 4). 
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Recalling now 1") we obtain 

± + 

37^ ^ X! k 

The function ^ , if it there exists, 
conclusion we lead if we consider 


-W< 1 C K> _ 
r L Z7* 

must satisfy eq. 16). To the same 


o 




As an application we can now restrict to the case C = C :e. 

is o 

We have C-(X) ^ and 

U.d.) =o 

1 di z k A 

This is a linear ordinary second, order equation. 

Its solution is of the form<y$ *6^and precisely 

+ xL (-kco -v L x ) 

4-k 




we obtain also 

«*•/<> ♦ y—z f l M *J 

In the previous analysis ^ is an arbitrary constant. If we let assume 
all the values 0 * 1 fa* 4 0 O we obtain ^as and (^complete sets of 

functions, which allow us to expande by integral whatever function. Each 
dynamic fenomenum in, a structure where eq.s 1") are valid and CT ^ 
can be analysed as, (see 17) and 17*)), 


If we consider 


<k 


4 = 


oo 



+ J % z (k) &*'!*+ 
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Obviously 18) is not the only way by which to expand by integral such a 
dynamic fenomenum, but this way allows us to consider component functions 
when the propagation speed at any x is well know. Such speed, as we know, 
is a fundamental datum in order to evaluate the speed of a reflected 
wave, in particular in the case of time depending restraint conditions. 
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SUMMARY OF THE SEPTEMBER 16 
TETHER DYNAMICS SIMULATION WORKSHOP 

Charles C. Rupp 

George C. Marshall Space Flight Center 

The objective of the Tether Dynamics Simu- 
lation Workshop is to provide a forum to 
discuss the structure and status of existing 
computer programs which are used to simulate 
the dynamics of a variety of tether appli- 
cations. The reasons for having different 
simulation models and how the simulations are 
verified is discussed. Recommendations made 
during a panel discussion regarding the dir- 
ection of future work are presented. 


INTRODUCTION 


This summary was prepared by a Tether Dynamics Simulation 
Review Committee consisting of the following members: 


Peter M. Bainum 

William A. Baracat 
Silvio Bergamaschi 
Paul A. Penzo 
Charles C. Rupp 


American As tronau t ical Soc iety /Howard 
Un ivers i ty 

General Research Corporation 

University of Padova, Italy 

Jet Propulsion Laboratory 

George C. Marshall Space Flight Center 


The workshop featured presentations on eleven simulation 
models and special presentations on the validation of the 
Tethered Satellite System Engineering Simulation at the 
Johnson Space Center and the dynamics flight experiments to 
be conducted on the first flight of the Tethered Satellite 
System. A panel composed of preselected workshop participants 
discussed five issues regarding future simulation validation 
activities. Approximately 120 people participated in the 
workshop. 


RECOMMENDED FUTURE SIMULATION DEVELOPMENT 


Comments that were received during the presentations of the 
simulation models are summarized as follows: 


1. Incorporate sensor hardware and observer measurement 
dynamics . 
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2. Improve simulation of tether transient thermodynamic 
effects such as rapid heating and cooling during 
terminator crossing. 

3. Resolve different approaches to atmospheric drag 
models. 

4. More experimental information is required on tether 
material damping and thermoelastic effects. 

5. Simplify simulations where possible. 

PANEL DISCUSSION TOPICS 

Five topics were selected for the panel discussion with the 
intent of provoking opposing viewpoints. However, little 
disagreement was found among the panel members and the mem- 
bers of the audience who participated. These topics, the 
names of the panel members leading the discussion, and the 
conclusions follow. 

TOPIC 1: 

Howard Flanders Martin Marietta Denver Aerospace 

Henry Wolf Analytical Mechanics Associates 

Should a universal simulation program be developed? 

a. General consensus is no. 

b. Desirable in theory, but not practical at this 
t ime . 

c. A universal program will be difficult to verify, 
costly to maintain, not likely to be widely and 
readily accepted. 

d. There may not exist such a program — various 
applications may require their own "best" program. 

e. A library of special purpose application programs 
or subroutines might be useful. 

f. These programs should be mini or micro-computer 
usable, well documented, and user-friendly. 

g. An alternative to a universal program development 
would be better documentation and use of existing 
programs. 

TOPIC 2: 

David A. Arnold Smithsonian As trophy s ical Observatory 

Silvio Bergamaschi University of Padova, Italy 

What accuracy is required in tether dynamics simulations? 

a. The simulation model should be selected based on 
the mission goals and the mission phase. 
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b. The influence of various factors should be studied 
individually to determine the magnitude of their 
effect. 

c. Effects should be included in the simulation 
according to their influence on the particular 
problem. 

d. Accuracy verification of the selected simulation 
should be made with a sample case of the complete 
simulation, and comparison of results with an 
independent program. 

e. Results depend not only on the accuracy of the 
computer algorithms, but on the accuracy of the 
physical parameters . 

TOPIC 3: 

Joe Carroll Energy Science Laboratory 

Charles C. Rupp George C. Marshall Space Flight Center 

What methods should be used to verify simulations? 

a. Mathematical confirmation can be used to compute 
simple cases whose results are known, compare with 
other simulation programs, and provide a self-test 
on the level of accuracy. 

b. Physical confirmation tests can be performed in a 
laboratory to determine tether physical character- 
istics, perform dynamics and control system tests, 
and test tether motion in a large vacuum chamber. 
The KC-135 can provide zero-G tests of tethers and 
system components. 

TOPIC 4: 

Enrico Lorenzini Smithsonian Astrophys ical Observatory 

Paul A. Penzo Jet Propulsion Laboratory 

Who should be responsible for simulation verification? 

a. The sponsoring agency should ultimately be 
responsible for verification. 

b. Possible options for the verification effort 
include in-house, unbiased third party, and/or 
a formal professional committee (e.g. AIAA-AAS 
connected ) . 

c. The methodology should include providing 
sufficient updated documentation for each model 
and establishing, case by case, a set of test 
runs . 
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TOPIC 5: 

Peter M. Bainum American As tronaut ical Soc iety/Howard 

University 

Bill Djinis NASA Headquarters 

Should a peer process be established? How? 

a. Numerous simulation programs have been developed 
under various contracts. 

b. The complexity and accuracy required of simulat- 
ions will depend on use (preliminary design, 
particular phase of mission support, experimental 
flight dynamics verification, future applicat- 
ions ) . 

c. Small new projects may not be able to afford 
development of their own simulations and could 
benefit from a group of expert peers. 

d. The panel recommends formation of a working group 
to study peer selection, availability, and 
financial support . 

RECOMMENDATIONS 

The final recommendations of the workshop are as follows: 

1. Establish a joint tether applications working group 
with the following objectives: 

a. Serve as a focal point for dynamics studies for 
tether applications . 

b. Advocate simulation validation. 

c. Maintain a simulation capabilities catalog. 

d. Perform periodic reviews of dynamics issues. 

e. Provide recommendations to management on resource 
requirements . 

2. Pursue ground and flight dynamics experimentation. 
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TETHER DYNAMICS SIMULATION WORKSHOP 
Tuesday, 16 September 1986 
PROGRAM 


7:45 
8: 15 

8:30 

9:00 

9:30 

10:15 

10:45 

11:20 

11:40 

12:00 

12:20 


Registration/Continental Breakfast 

Introduction 
Chris Rupp 

NASA Marshall Space Flight Center 

William A. Baracat 
General Research Corporation 


Simulation Descriptions 

Tether Dynamics Simulation 
Vinod J. Modi 

University of British Columbia, Canada 

Arun K. Misra 

McGill University, Canada 

HIFITS Simulation 
Howard Flanders 

Martin Marietta Denver Aerospace 

Tethered Satellite System on the Systems Engineering Simulator 
Ron Humble 

Lockheed Missiles and Space Company 
David Harshman 

Rockwell Space Operations Company 
GTOSS 

David D« Lang 

David D. Lang Associates 

Tethered Satellites Simulation 

John R. Glaese 

Control Dynamics Company 

Skyhook Program 
David A. Arnold 

Smithsonian Astrophysical Observatory 
Slack3 

Gordon E. Gullahorn 

Smithsonian Astrophysical Observatory 

Artificial Gravity Laboratory 
Enrico Lorenzini 

Smithsonian Astrophysical Observatory 


Lunch/Computer Simulation Demonstrations 
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PM 

1:30 


Simulation Descriptions - Continued 


2:00 


2:45 


3:15 


3:30 


4:00 


Elasticity Effects in Tether Dynamics 
Silvio Bergamaschi 
University of Padova, Italy 

TSS Subsatellite Attitude Dynamics and Control Laws 
Verification Programs 
Floriano Venditti 

Aeritalia Space Systems Group, Italy 

Validation of TSS Engineering Simulation 
Zachary Galaboff 

NASA Marshall Space Flight Center 

TSS-1 Dynamics Flight Experiments 
Gordon E. Gullahorn 

Smithsonian Astrophysical Observatory 

Tether Simulation Design for Mission Planning and Analysis 

Richard Deppisch 

NASA Johnson Space Center 


Yashvant Jani 
LinCom Corporation 

Panel Discussion on Future Validation Activities 


Topic 1: Howard Flanders, Henry Wolf: Should a single 

universal simulation program on tether dynamics be 
built? 


Topic 2: Dave Arnold, Silvio Bergamaschi: How accurately must 

simulations be able to compute tether system 
dynamics? 

Topic 3: Joe Carroll, Chris Rupp: What verification methods 

should be employed? 


Topic 4: Enrico Lorinzini, Paul Penzo: Who should be 

responsible for simulation verification? 

Topic 5: Peter Bainum, Bill Djinis: Should a peer process be 

established? 


5:30 Adjournment 

6:00 Executive Committee Meeting - Report to be presented Friday 

10 : 00 . 


340 



APPENDIX B 


WORKSHOP ATTENDEES 


341 


Randall S. Aadland 
NASA Lewis Research Center 
NASA— CASE Co-op Student 
21000 Brookpark Road 
Cleveland, OH 44135 
(216) 433-6134 

Adrian Abraham 

Bethesda-Chevy Chase Space Experiment Club 

Club Leader 

6406 Stratford Road 

Chevy Chase, MD 20815 

(301) 654-0127 

Dale T, Adams 
Texas A&M University 
Military Aerospace Association 
PO Box 6500 

College Station, TX 77844 
(409) 260-7140 

Charles Alexander 
Energy Science Labs 
PO Box 85608 

San Diego, CA 92138-5608 
(619) 455-4688 

Norman Alexander 

General Electric Space Div. 

P0 Box 8555 Bldg. 11 
Philadelphia, PA 19101 
(215) 962-2123 

Kyle T. Alfriend 
General Research Corporation 
7655 Old Springhouse Road 
McLean, VA 22102 
(703) 893-5900 

Marc Allen 

Computer Technology Associates 
11847 Canon Blvd., Suite 3 
Newport News, VA 23606 
(804) 873-2021 

Sean Allen 
Field Box 183 
Wapello, Iowa 52653 


343 


msgeding page blank not film® 


John L. Anderson 
NASA Headquarters 

Office of Aeronautics and Space Technology 
Mail Code RS 
Washington, DC 20546 
(202) 453-2756 

Ettore Antona 

Politecnico Di Torino 

Dipartimento Ingegneria Aerospaziale 

Corso Duca degli Abruzzi, 24 

10146 Torino, Italy 

39-11-5566807 

David A. Arnold 

Smithsonian Astrophysical Observatory 

60 Garden Street 

Mail Stop 59 

Cambridge, MA 02138 

(617) 495-7269 

Wayne L. Bailey 
Teledyne Brown Engineering 
300 Sparkman Drive 
MS 52 

Huntsville, AL 35807 
(205) 532-2829 

Peter M* Bainum 
Howard University 

Department of Mechanical Engineering 
Washington, DC 20059 
(202) 636-6612 

William Baracat 

General Research Corporation 

Aerospace Group 

7655 Old Springhouse Road 

McLean, VA 22102 

(703) 893-5900 

Ivan Bekey 
NASA Headquarters 
Director, Advanced Programs 
Code MT 

Washington, DC 20546 
(202) 453-1147 

Silvio Bergamaschi 
University of Padua 
Institute of Applied Mechanics 
Via Venezia 1 
35131 Padua, Italy 
39-49-8071033 


344 



Gianfranco Bevilacqua 
Aeritalia 

Space Systems Group 
Corso Marche 41 
10146 Torino, Italy 
39-11-725089 

Giannandrea Bianchini 
University of Padua 
Via Venezia 1 
35131 Padua, Italy 
39-49-20198 

Karen D. Blumentritt 
Rockwell International 
1840 NASA Road One 
Houston, TX 77058 
(713) 483-2627 

Andrew E. Boeck 
6558 Idlewild Road 
Apartment E 
Charlotte, NC 28212 
(504) 242-0714 

Alexander Boschitsch 

Princeton University 

Dept, of Mechanical Engineering 

D-209 Eng. Quadrangle 

Princeton, NJ 08544 

(609) 452-5674 

Edward J. Brazill 
NASA Headquarters 
Code MTS 

Washington, DC 20546 
(202) 453-1155 

Cyrus L. Butner 

General Research Corporation 

Aerospace Group 

7655 Old Springhouse Road 

McLean, VA 22102 

(703) 893-5900 

Michael A. Calabrese 
NASA Headquarters 
Code EED 

Washington, DC 20546 
(202) 453-1676 


345 


Nellie N. Carr 
NASA Johnson Space Center 
AST Flight Systems Operations 
Code DH64 

Houston, TX 77058 
(713) 483-4966 

Joseph A. Carroll 
Energy Science Labs 
P0 Box 85608 

San Diego, CA 92138-5608 
(619) 455-4688 

Brian F. Class 
Fairchild Space Company 
Senior Engineer, ACS Section 
2031 Century Blvd. 

Germantown, MD 20874 
(301) 428-6245 

Mario Cosmo 

Smithsonian Astrophysical Observatory 
60 Garden Street 
Cambridge, MA 02138 
(617) 495-7433 

Capt. Carl E. Crockett 
United States Air Force 
Air Force Space Command/DOA 
HQ Space Command 
Peterson AFB 

Colorado Springs, CO 80914 
(303) 554-2851 

Eric Dahlstrom 
7630 Tomlinson Avenue 
No. 35 

Cabin John, MD 20818 
(301) 229-6599 

Richard M. Depplsch 
NASA Johnson Space Center 
Mission Planning and Analysis 
Houston, TX 77058 
(713) 483-4457 

Lamont DiBiasi 

Fairchild Space Company 

Director, NASA Business Development 

20301 Century Boulevard 

Germantown, MD 20874-1181 

(301) 428-6363 


346 



William Djinis 
NASA Headquarters 
Code MTF 

Washington, DC 20546 
(202) 453-1157 

Charles R. Eastwood 
NASA Headquarters 
Code EED 

Washington, DC 20546 
(202) 453-1723 

Pamela Edwards 
AIAA 

1633 Broadway 
New York, NY 10019 
(212) 408-9778 

Marinella Ercoli 
PSN/CNR 

Viale Regina Margherita 202 
00198 Roma, Italy 
39-6-8445107 

Ruying Fan 

Howard University 

Dept, of Mechanical Engineering 

Washington, DC 20059 

(202) 636-7124 

Joseph V. Fedor 

NASA Goddard Space Flight Center 
Code 712 

Greenbelt , MD 20771 
(301) 286-2088 

Mr. Festa 
Aeritalia 

Space Systems Group 
Corso Marche 41 
10146 Torino 
Italy 

39-11-333271 
Howard Flanders 

Martin Marietta Denver Aerospace 
Mail Stop S8071 
PO Box 179 
Denver, CO 80201 
(303) 977-5951 


347 


Kevin Flannery 
General Electric RSL 
3129 Chestnut St* 

Room 5450M 

Philadelphia, PA 19101 
(215) 823-3759 

Thomas W. Flatley 

NASA Goddard Space Flight Center 

Code 712 

Greenbelt , MD 20771 
(301) 286-5336 

Zachary J. Galaboff 

NASA Marshall Space Flight Center 

Code ED 13 

Marshall Space Flight Center, AL 
(205) 544-1446 

Marco Galantino 
PSN/CNR 

Viale Regina Margherita 202 
00198 Roma, Italy 
39-6-4767249 

William H* Ganoe 
University of Arizona 
S.I*E. Department 
Tucson, AZ 85721 
(602) 621-6562 

James Gearhart 

777 West Middlefield Road 

#134 

Mountain View, CA 94043 
(408) 742-8085 

Mireille Gerard 
AIAA 

Administrator, Corporate and 
International Programs 
1633 Broadway 
New York, NY 10019 
(212) 408-9770 

Francesco Giani 
Aeritalia 

Space Systems Group 
Corso Marche 41 
10146 Torino, Italy 
39-11-3332235/721345 


35812 


348 



John R. Glaese 
Control Dynamics Company 
Office Park South 
Suite 304 

600 Boulevard South 
Huntsville, AL 35802 
(205) 882-2650 

Ron Glickman 
Ball Aerospace 
P.0. Box 1062 
Boulder, CO 80306 
(303) 939-4275 

Filippo Graziani 
Universita Di Roma 
Associate Professor 
Via Lariana 5 
Roma 00199 
Italy 

39-6-858333 

Luciano Guerriero 
PSN/CNR 

Viale Regina Margherita 202 
00198 Roma, Italy 
39-6-4767249 

Gordon E* Gullahorn 

Smithsonian Astrophysical Observatory 

Mail Stop 59 

60 Garden Street 

Cambridge, MA 02138 

(617) 495-7419 

Kevin Hames 

University of Michigan 
5594 Northcote 
West Bloomfield, MI 48011 
(313) 626-0509 

Calvin D. Harr 

Martin Marietta Denver Aerospace 
7775 S. Trenton Ct. 

Englewood, CO 80112 
(303) 977-5451 

James K. Harrison 

NASA Marshall Space Flight Center 

Mail Code PS04 

Marshall Space Flight Center, AL 35812 
(205) 544-0629 


349 


Dave Harshman 

Rockwell Space Operations Co. 

600 Gemini Avenue 
Houston, TX 77058 

(713) 282-2789 

William S. Hayes 
McDonnell Douglas Astronautics 
14651 Mimosa Lane 
Tustin, CA 92680 

(714) 896-1338 

Xiaohua He 
Stanford University 
Aeronautics/Astronautics Dept. 

Stanford, CA 94305 
(415) 725-3297 

Gail Heim 

System Planning Corporation 
1500 Wilson Blvd. 

Suite 1202 

Arlington, VA 22209 
(703) 841-2880 

Ronald W. Humble 

Lockheed Engineering & Management Services 

Mail Code B12 

2400 NASA Road 1 

Houston, TX " 77058 

(713) 333-6285 

Paul Ibanez 

ANCO Engineers, Inc. 

9937 Jefferson Blvd. 

Culver City, CA 90232-3591 
(213) 204-5050 

Rodica Ionasescu 
Jet Propulsion Laboratory 
4800 Oak Grove Drive 
Pasadena, CA 91109 
(818) 354-0445 

Yashvant Jani 
LinCom Corporation 
18100 Upper Bay Road 
Suite 208 

Houston, TX 77058 


350 



Philip 0. Jarvinen 

Massachusetts Institute of Technology 
244 Wood Street, Rm. D-382 
Lexington, MA 02173 
(617) 863-5500 

Francis J. Kelly 

Naval Research Laboratory 

Section Head, Long Wave Radio 

4555 Overlook Avenue 

Washington, DC 20375 

(202) 767-2941/(202) 767-2891 

Joseph C. Kolecki 
NASA Headquarters 
Code MTS 

Washington, DC 20546 
(202) 453-1881 

Meredith C. Kurz 
4101 Longfellow Street 
Hyattsville, MD 20781 
(301) 277-5682 

David D. Lang 
David D. Lang Associates 
2222 70th Avenue, SE 
Mercer Island, WA 98040 
(206) 236-2579 

James R. Lease 
NASA Headquarters 
Code MSD 

Washington, DC 20546 
(202) 453-1576 

John F. L. Lee 
Honeywell Inc. 

Mail Station 809-5 
13350 U.S. Highway 19, South 
Clearwater, FL 33546 
(813) 539-5667 

Carolina Leonard 
Adroit Systems 
809 N. Royal Street 
Alexandria, VA 22314 
(703) 684-2900 

Sue Leszkiewicz 

NASA Goddard Space Flight Center 
Code 554.1 

Greenbelt , MD 20771 
(301) 286-3155 


351 


Alejandro Levi 
ANCO Engineers, Inc. 

9937 Jefferson Blvd. 

Culver City, CA 90232-3591 
(213) 204-5050 

John Lewak 
Bristol Aerospace 
P0 Box 874 
Winnipeg R3C 2S4 
Canada 

(204) 775-8331 

Lin Liangdong 
Howard University 
2225 N Street, NW 
Apartment 404 
Washington, DC 20037 
(202) 636-7124 

Enrico C. Lorenzini 

Smithsonian Astrophysical Observatory 

Mail Stop 59 

60 Garden Street 

Cambridge, MA 02138 

(617) 495-7211 

Andrea Lorenzoni 
PSN/CNR 

Viale Regina Margherita 202 
00198 Roma, Italy 
39-6-4767259 

Charles A. Lundquist 

University of Alabama in Huntsville 

Research Institute 

Director of Research/CMDS 

Room M-65 

Huntsville, AL 35899 

(205) 895-6620 

Gianfranco Manarini 
PSN/CNR 

Viale Regina Margherita 202 
00198 Roma, Italy 
39-6-4767249 

Chris A. Markland 

European Space Agency 

2200 AG 

Noordwijk 

The Netherlands 

31-17198655 


352 



Pietro Merlina 
Aeritalia 

Space Systems Group 
Corso Marche 41 
10146 Torino 
Italy 

39-11-3332356 
Paul A. Miller 

RCA Astro-Electronics Division 
PO Box 800 

Princeton, NJ 08543 
(609) 426-3857 

William Miller 
Aeritalia 

Space Systems Group 
Corso Marche 41 
10146 Torino, Italy 
39-11-725089 

Arun K. Misra 
McGill University 

Department of Mechanical Engineering 
817 Sherbrooke Street, W. 

Montreal, Quebec, Canada H3A 2K6 
(514) 392-5480 

Antonio Moccia 
Istituto di Gasdinaraica 
Universita di Napoli 
Pie Tecchio 80 
08125 Napoli, Italy 
39-81-7682158 

Vinod J. Modi 

University of British Columbia 
Department of Mechanical Engineering 
2324 Main Mall 
Vancouver, B.C. 

Canada V6T 1W5 
(604) 228-2914 

Dean S. Monitor 

Martin Marietta Aerospace 

Michoud Operations 

Dept. 3514 

P.0. Box 29304 

New Orleans, LA 70189 

(504) 255-2092 


353 


Maria B. Montini 
PSN/CNR 

Viale Regina Margherita, 202 
00198 Roma, Italy 
39-6-4767261 

Claude Nicollier 

NASA Johnson Space Center 

ESA 

Office Code CB 
Houston, TX 77038 
(713) 483-4513 

John Oldson 

Energy Science Laboratories 
P0 Box 85608 

San Diego, CA 92138-5608 
(619) 455-4688 

Michael Patterson 
NASA Lewis Research Center 
Mail Stop 500-220 
21000 Brookpark Road 
Cleveland, OH 44135 
(216) 433-2405 

Jess A. Pedrena 
Academy of Aeronautics 
La Guardia Airport 
530 E. 20th Street 
New York, NY 10009 
(212) 254-6536 

Paul A* Penzo 
Jet Propulsion Laboratory 
Mail Code 301-170U 
4800 Oak Grove Drive 
Pasadena, CA 91109 
(818) 354-6162 

J* David Powell 
Stanford University 
Department of Aeronautics and 
Astronautics 
Durand Building 
Stanford, CA 94305 
(415) 723-3425 

Christopher R* Purvis 
Jet Propulsion Laboratory 
Mail Code 156-246 
4800 Oak Grove Drive 
Pasadena, CA 91109 
(818) 354-4986 


354 



Terrence G. Reese 
General Research Corporation 
Director, Aerospace Group 
7655 Old Springhouse Road 
McLean, VA 22102 
(703) 893-5900 

Charles C. Rupp 

NASA Marshall Space Flight Center 
Mail Code PS04 

Marshall Space Flight Center, AL 35812 
(205) 544-0627 

Cal Rybak 
Ball Aerospace 
P.0. Box 1062 
Boulder, CO 80306 
(303) 939-5462 

Jean Sabbagh 
PSN/CNR 

Viale Regina Margherita 202 
00198 Roma, Italy 
39-6-4767247 

Lester Sackett 

Charles Stark Draper Laboratory, Inc. 

Mail Stop 4-C 

555 Technology Square 

Cambridge, MA 02139 

(617) 258-2283 

Chikatoshi Satoh 
Nihon University 
7-24-1, Narashinodai 
Funbashi-city, Chiba-pref 
274 Japan 
81-47-466-1111 

Claire Savaglio 
University of Michigan 
Aerospace Engineering 
3217 Harvard 
Royal Oak, MI 48072 
(313) 288-5668 

Ralph Schinnerer 
Ford Aerospace 
10440 State Highway 83 
Colorado Springs, CO 80908 
(303) 594-1367 


355 


Silvano Sgubini 
Universita Di Roma 
Dipartimento Aerospaziale 
Via Salaria 851 
00100 Roma, Italy 
39-6-8120529 

Stanley D. Shawhan 
NASA Headquarters 

Branch Chief, Space Plasma Physics 
Mail Code EES 
Washington, DC 20546 

(202) 453-1676 

Anna Sinopoli 
Universita di Venice 
Assistant Professor 
S. Croce 191 
30125 Venice, Italy 
39-49-8071033 

Emil Soderberg 
Analysis & Technology 
258 Bank Stret 
New London, CT 06320 

(203) 444-0827 

Carl I. Soderland 
NASA Johnson Space Center 
Code EH24 
Houston, TX 77058 
(713) 483-4405 

Thomas H. Stengle 

NASA Goddard Space Flight Center 

Mail Code 554.1 

Greenbelt, MD 20771 

(301) 286-6354 

Stephen W. Strom 
TRW Defense Systems Group 
PO Box 58327 
Houston, TX 77258 
(713) 333-3133 

Thomas D. Stuart 
NASA Headquarters 
Office of Space Flight 
Mail Code MST 
Washington, DC 20546 
(202) 453-1562 


356 



Daniel Swanson 
Princeton University 
MAE Department 
Eng. Quadrangle 
Princeton, NJ 08544 
(609) 452-6083 

Frank R. Swenson 
Tri-State University 
Professor, Mechanical Engineering 
Angola, IN 46703-0307 
(219) 665-3141, ext. 256 

Giorgio Tacconi 

University of Genova 

Dipartimento Dibe Universita Genova 

Via dell'Opera Pia 11 A 

16145 Genova, Italy 

39-10-311811 

John M. Van Pelt 

Martin Marietta Denver Aerospace 
Mail Stop S8071 
PO Box 179 
Denver, CO 80201 
(303) 977-4888 

Floriano Venditti 
Aeritalia 

Space Systems Group 
Dynamic Systems Director 
Corso Marche 41 
10146 Torino, Italy 
39-11-3332752 

Sergio Vetrella 
Universita di Napoli 
Istituto Gasdinamica 
P.le Tecchio 80 
08125 Napoli, Italy 
39-81-7682159 

Andy Von Flotow 

Massachusetts Institute of Technology 
Department of Aeronautics and Astronautics 
7755 Massachusetts Avenue 
Cambridge, MA 02139 
(617) 253-4865 


357 


Roger Wacker 

McDonnell Douglas Astronautics Co. 

Mail Code T5K 

16055 Space Center Blvd. 

Houston, TX 77062 
(713) 280-1772 

Justin Walls 

University of Alabama at Huntsville 
Huntsville, AL 35899 
(205) 895-6126 

Douglass Wheeluck 

University of Alabama at Huntsville 
Huntsville, AL 35899 
(205) 895-6126 

Roger Williamson 
Stanford University 
STARLab/SERA 
Durand Bldg., Room 331 
Stanford, CA 94305-4055 
(415) 723-9774 

Henry Wolf 

Analytical Mechanics Associates, Inc. 
17 Research Drive 
Hampton, VA 23666-1398 
(804) 865-0944 

Holger Wolff 

Messerschmitt Boelkow Blohn GMBH 
Postfach 801169 
8012 Ottobrunn 
Muenchen, West Germany 
49-84-60007762 

William Woodis 

Martin Marietta Denver Aerospace 
Mail Stop S-8073 
PO Box 179 
Denver, CO 80201 
(303) 977-6452 

Italo Zonin 

Embry-Riddle Aeronautical University 
327 Wilder Blvd. 

Unit 6-104 

Daytona Beach, FL 32014 
(904) 252-4976 


358 



NASA 


2. Government Accession No. 


3. Recipient's Catalog No. 


4. Title and Subtitle 

Tether Dynamics Simulation 


5. Report Date 
February 1987 

6. Performing Organization Code 
MTS 


7. Author(s) 


8. Performing Organization Report No. 


10. Work Unit No. 


9. Performing Organization Name and Address 

NASA Office of Space Flight, Advanced Programs 
Washington, DC 20546 

and 

George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 


11. Contract or Grant No. 


13. Type of Report and Period Covered 
Conference Publication 

14. Sponsoring Agency Code 


15. Supplementary Notes 

The proceedings of the International Conference on Tethers in Space, 
published under the auspices of the American Astronaut ica 1 Society under 
the title Volume 62, Advances in the As t ronau t ica 1 Sciences; Tethers in 
Space , are available from Univelt, Inc., P.0. Box 28130, San Diego, CA 92128. 

16. Abstract 

The proceedings of the first Tether Dynamics Simulation Workshop are 
presented. The workshop was held at the Hyatt Regency Crystal City Hotel 
in Arlington, Virginia, September 16, 1986, in conjunction with the 
International Conference on Tethers in Space. The objective of this work- 
shop was to provide a forum for the discussion of the structure and status 
of existing computer programs which are used to simulate the dynamics of 
a variety of tether applications in space. A major topic of the workshop 
was the purpose of having different simulation models and the process of 
validating them. Guidance on future work in these areas was obtained 
from a panel discussion; the panel was composed of resource and technical 
managers and dynamic analysts in the tether field. The conclusions of 
this panel are also presented in this document. 


17. Key Words (Suggested by Author(s)) e 1 ec t ro dynamic 18. Distribution Statement 
tethers space shuttle tethers 

dynamics space station Unclassified - Unlimited 

computer simulation 

tethered satellite system c , . _ r* 77 

, , 7 Subject Category 37 

tether dynamics 

19. Security Classif. (of this report) 20. Security Classif. (of this page) 21. No. of pages 

Unclassified Unclassified 368 


22. Price 


NASA-Langley, 1987 



